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Chapter 1 

Introduction 


Flywheels offer great promise for reducing the mass and extending the life of spacecraft; they store 
more energy per unit of mass and last significantly longer than chemical batteries. Moreover, 
flywheels can simultaneously store energy and exert torque on a spacecraft, making it possible 
for one system of flywheels to replace two separate systems typically used for energy storage and 
attitude control. When the mass of the two conventional systems is taken into account, the specific 
energy of flywheel systems is expected to be 5 to 10 times greater, according to Ref. [1]. The 
attitude control system typically represents 11% of the mass of a spacecraft, and batteries make up 
6% of the mass; replacing 17% of a spacecraft’s mass with a flywheel system whose mass is 1.7% 
would lead to a 15% reduction in total spacecraft mass. Secondary benefits occur as well; since 
flywheels have higher system level efficiencies than batteries, a reduction in solar array size and mass 
becomes possible, and reboost propellant can be reduced because the smaller arrays produce less 
drag. Flywheel systems are expected to last 15 years (Ref. [1]) or more whereas typical batteries 
last only 5 years. The greatest advantage of flywheels over batteries accrues in low earth orbit 
where eclipse happens more frequently and for larger fractions of an orbit than in higher orbits; 
repeated charge and discharge cycles, and high depth of discharge significantly degrade batteries 
over time. 

Flywheel energy storage technology can be employed in the transportation, utility, and man- 
ufacturing industries; potential NASA applications include spacecraft, launch vehicles, aircraft 
power systems, uninterruptable power supplies, and planetary rovers. Flywheels are recognized as 
a critical space technology in NASA’s recent Advanced Cross-Enterprise Technology Development 
research announcement, Ref. [2], which called for “high-risk, high-payoff technology advances.” The 
Advanced Power and On-Board Propulsion section of the announcement requested proposals for 
“energy storage . . . techniques that dramatically reduce spacecraft bus mass, volume, and cost”; 
specifically, “flywheel energy storage system component improvements . . . that enhance their . . . 
integration with spacecraft attitude control functions.” To date, flywheel energy storage devices 
have not yet been flown onboard an actual spacecraft. 

The Agency’s interest in flywheel technology spans several decades, including two workshops 
devoted to flywheels in 1983 and 1984, two workshops dealing with magnetic suspension held at 
Langley Research Center (LaRC) in 1988 and 1990, and flywheel workshops in 1997, 1998, 1999, 
2000, and 2002 co-sponsored by Glenn Research Center (GRC), where an active development 
program exists. Under contracts administered by LaRC, Integrated Power and Attitude Control 
Systems (IPACS) were evaluated in the mid 1970’s as reported in Ref. [3], and considered for the 
Space Station in the mid 1980’s as summarized in Ref. [4]. 
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1.1 International Space Station Flywheels 


For several years the International Space Station (ISS) Payloads Office funded a program at GRC 
to develop flywheel energy storage technology. The project would have added significantly to the 
capacity for energy storage onboard the Station, and reduced or eliminated the cost and time 
required to replace chemical batteries. 

Each device in the ISS Flywheel Energy Storage System (FESS) was to consist of two counter- 
rotating rotors placed in vacuum housings, and levitated with magnetic bearings. Motor-generators 
would connect the rotors to the existing electrical power system so that they could store energy 
when available from the photovoltaic arrays, and supply energy when needed. Each rotor, made up 
of a metallic hub and a rim of composite material, would spin at angular speeds ranging between 
18,000 and 60,000 rev/min. The FESS could have been made up of as many as 48 flywheel pairs, 
replacing all of the batteries on the Assembly Complete configuration of ISS. Unfortunately, plans 
called for FESS to be used only for energy storage, not for attitude control. 

The attitude of ISS is controlled with 4 Control Moment Gyroscopes (CMGs), together with a 
Reaction Control System (RCS). The CMGs were activated following installation of the laboratory 
module Destiny in February, 2001, and 1 CMG experienced a bearing failure in June, 2002, after 
less than 16% of the expected 8^ year performance life. CMGs are able to control attitude without 
requiring expenditure of RCS propellant most of the time; however, in certain situations such as 
Orbiter rendezvous and docking, the CMGs will become saturated. Desaturation with the RCS 
consumes propellant, which must then be resupplied from Earth at a cost of $10,000 to $20,000 per 
pound. 

Forty-eight pairs of FESS flywheels would have possessed a total useable angular momentum 
(Ref. [5]) of approximately 120,000 N-rn-s; only 50% to 75% of this figure would have been required 
for energy storage, leaving 30,250 to 60,500 N-m-s available for attitude control, or 1.5 to 3 times 
the amount available from the current set of 4 CMGs. The FESS, or a similar system, could thus 
provide a backup system for momentum exchange that could serve even in the event of failure of all 
CMGs. Used to assist the CMGs, a flywheel system could allow attitude maneuvers that are faster 
or of larger excursions than would be possible with the CMGs alone, and better accommodation of 
activities that require CMG momentum, such as Remote Manipulator System operations with large 
payloads, centrifuge spin-up and spin-down, and the Alpha Magnetic Spectrometer experiment with 
its strong permanent magnets. Steady-state attitude motion could be reduced, thus improving the 
quasi-steady microgravity environment. It may even be possible to reduce or eliminate expenditure 
of propellant required occasionally to desaturate the CMGs. 

1.2 Existing Control Laws 

A numerical investigation of the merits of using the FESS to assist the CMGs requires a feedback 
control law designed for CMGs and flywheels used together as an integrated set of effectors. The 
current CMG control law was originally developed by Wie et al. in Ref. [6] and refined by Harduvel 
as described in Ref. [7]; it minimizes a cost function involving spacecraft attitude and angular 
speed, and CMG angular momentum. A CMG steering law, such as the one developed by Kennel 
in Ref. [8], determines the speeds of the two gimbals (in which each constant speed CMG rotor is 
mounted) needed to produce the torque requested by the control law. As a natural extension of 
the present approach, we seek a new control law derived from a cost function that includes flywheel 
angular momentum in addition to the aforementioned quantities. Also needed is a “flywheel steering 
law,” a counterpart to the CMG steering law that will determine the motor-generator torque to be 
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applied to each member of a flywheel pair such that rotational kinetic energy is stored or discharged 
in the required manner, and the net torque requested by the control law is produced simultaneously. 

A review of the literature does not reveal any existing three-axis control laws for earth-pointing 
spacecraft using flywheels and CMGs together, or even flywheels alone, where attitude control, 
momentum management, and power management are addressed in a unified way. As mentioned 
earlier, the current ISS attitude control system is designed only for use of CMGs (and CMG 
desaturation with RCS); no provisions exist for using flywheels or other momentum exchange 
devices. In Ref. [3], Notti, Cormack, and Klein give a sketch of a control law and an energy 
distribution law; however, this work is not applicable primarily because each flywheel rotor is 
assumed to be supported by two gimbals and the ISS FESS did not contain any gimbals. In 
addition, the control law lacks flywheel angular momentum as a feedback parameter. Recent work 
on control laws for integrated power and attitude control systems deals either with sets of four or 
more flywheels whose spin axes are non-collinear (Refs. [9]— [11] ) , or with sets of variable-speed, 
single-gimbal control moment gyroscopes (Refs. [12]— [14] ) ; neither of these types of configurations 
are directly applicable to the counter-rotating flywheel pair arrangement of the FESS. Varatharajoo 
and Fasoulas develop control laws in Ref. [15] for a spacecraft using a counter-rotating flywheel 
configuration; however, they only consider a pitch-axis controller. 

Hall’s control law, proposed in Ref. [9] , is an open-loop scheme (rather than a feedback law) for 
performing large-angle attitude maneuvers. It does not account for gravitational and aerodynamic 
torques which have a significant effect on the motion of ISS, and therefore can not be used for 
maintaining torque equilibrium attitude, the primary job of the CMGs. Hall introduces a flywheel 
steering law based on a matrix pseudo-inverse; it is applied in each of Refs. [10]— [14] . Tsiotras, 
Shen, and Hall employ Lyapunov stability theory in Ref. [10] to develop a feedback control law 
which performs well in simulations involving disturbance torques; however, flywheel momentum 
is managed by expenditure of propellant. Costic et al. develop in Ref. [11] a nonlinear controller 
which includes an adaptive scheme for estimating the mass distribution of the spacecraft, but they 
do not address momentum management. 

In Ref. [12], Fausz and Richie extend the work of Hall to a nonlinear feedback controller appli- 
cable to a set of variable-speed, single-gimbal control moment gyroscopes. Together with Tsiotras, 
they continue their discussion in Ref. [13] and present simulation results, but momentum manage- 
ment is not addressed in either of the two papers. Yoon and Tsiotras develop an adaptive nonlinear 
control law in Ref. [14], and incorporate wheel-speed equalization to reduce the possibility of sin- 
gularities and keep the wheel speeds within acceptable limits. Numerical simulation results show 
that attitude and power profiles can be tracked even when the spacecraft inertia properties are 
unknown. Of all the papers mentioned heretofore, Ref. [14] is the only one in which attitude con- 
trol, momentum management, and power management for flywheels is considered in an integrated 
fashion. 

It is important to note that Refs. [9] — [15] fail to take into account damping torque exerted by 
the spacecraft and a flywheel rotor on each other; in practice, this will cause the actual rotational 
kinetic energy possessed by the flywheels to differ from the required amount. None of these works 
include any feedback of errors in flywheel power or kinetic energy, something which must be done 
under realistic conditions. 

1.3 Overview 

As a prerequisite for the design of laws for controlling attitude and managing angular momentum 
of flywheels and CMGs, one must have in hand equations of motion for a spacecraft containing 
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these devices. In Chapter 2 of the material that follows, a complete, nonlinear set of equations 
of motion for such a system is derived in vector-dyadic form, and subsequently compared to some 
others existing in the literature. Two configurations of special interest are then examined. In 
the absence of CMGs, the system forms a complex gyrostat; scalar equations of motion are given 
for the case of three orthogonally mounted pairs of flywheels where the spin axes of the members 
of each pair are parallel. Another special case occurs when flywheels are absent; under certain 
assumptions, the equations are shown to correspond to those used by Wie et al. in Ref. [6] for a 
spacecraft with CMGs. These assumptions made in connection with the CMGs are then applied 
to the full system, resulting in approximate equations which are subsequently linearized about an 
Earth-pointing motion, and nondimensionalized for use in control law design. 

Two flywheel “steering laws” are developed in Chapter 3 such that torque commanded by an 
attitude control law is achieved while energy is stored or discharged at the required rate. The first 
law is based upon a suggestion by Hall in Ref. [9]. It requires the formation of a matrix pseudo- 
inverse; the total power requirement is provided by the entire flywheel system, with no assumptions 
regarding the distribution of the requirement among the three flywheel pairs. The second steering 
law follows from the assumption that the power requirement is divided evenly among the three 
flywheel pairs, and does not require construction of a pseudo-inverse. 

In Chapter 4 designs are given for control laws, obtained with the Linear Quadratic Regulator 
technique, that allow CMGs and flywheels to be used as an integrated set of effectors for maintaining 
torque equilibrium attitude and managing momentum, or for holding a specified attitude until 
momentum saturation occurs. The control algorithms compensate for damping torque, exerted by 
the flywheel bearings, through feedback of error in rotational kinetic energy, that is, the difference 
between actual energy of the flywheels and the energy they should have based on the required 
power. In numerical simulations of nonlinear motion, several variations of the linear control laws 
are seen to perform well in achieving torque equilibrium attitude when initial conditions of the 
spacecraft differ from the steady-state conditions, and the spacecraft is subject to gravity-gradient 
torque together with a prescribed time history of aerodynamic torque. Angular momentum of 
flywheels and CMGs is minimized, flywheels store or discharge energy as commanded, and flywheel 
angular speeds are kept within appropriate limits. 

Instead of continuously varying the spacecraft attitude to manage angular momentum, the 
controllers can be requested to keep the spacecraft in a specified orientation until the angular 
momentum capacity of the actuators is exceeded, at which time thrusters can be used for the 
remainder of the operation. Hence, the controllers designed for momentum exchange devices have 
been supplemented with a minimum-fuel jet selection scheme that enables thrusters to provide the 
required torque. Flywheels can be used to aid the CMGs in holding attitude, resulting in a delay 
for the time of CMG momentum saturation and onset of thruster firings, and thus in a reduction of 
the propellant expended. Simulation results indicate a certain amount of propellant can be saved 
during operations leading up to docking between ISS and another spacecraft. 

Nonlinear control algorithms for performing large-angle reorientations or slew maneuvers are 
presented in Chapter 5. These control algorithms are based on feedback linearization and Lyapunov 
stability theory, and they take into account constraints such as slew rate limits and maximum 
control torques. Good performance is observed in simulations; however, the momentum capacity of 
ISS CMGs, together with the expected capacity of ISS flywheels, would be insufficient to perform 
most Space Station maneuvers without expending propellant. 

Chapter 6 concludes the report with a summary of the results, as well as some suggestions for 
future research. 
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Chapter 2 

Equations of Motion for Spacecraft 
with Flywheels and CMGs 


One must have in hand equations governing the motion of a spacecraft possessing flywheels and 
CMGs if laws for controlling the motion are to be crafted, and subsequently exercised in numerical 
simulations. Kane’s method is used in Sec. 2.1 to derive a set of nonlinear dynamical equations in 
vector-dyadic form; they are general with regard to the number and orientation of flywheel rotors, 
and the number of CMG rotors and gimbals. Certain terms of the equations are shown to be in 
agreement with work previously published by others. In Sec. 2.2 a set of twelve scalar equations is 
obtained by applying the generic relationships to the special case of a complex gyrostat consisting 
of a base body and three pairs of flywheels mounted in orthogonal directions. When only one 
flywheel is present the system becomes a simple gyrostat, and the resulting equations for rotational 
motion are shown to correspond to expressions given by Kane, Likins, and Levinson in Ref. [16]. A 
second configuration of particular interest arises when CMGs are the only spacecraft appendages; 
the equations used by Wie et al. in Ref. [6] are shown in Sec. 2.3 to follow from the general 
ones given in Sec. 2.1, along with several reasonable, clearly stated assumptions. The exact and 
approximate relationships presented in Secs. 2.2 and 2.3 respectively are combined in Sec. 2.4 to 
form approximate equations for a spacecraft with flywheels and CMGs, which are linearized in 
Sec. 2.5 and nondimensionalized in Sec. 2.6 in preparation for design of linear control laws. Finally, 
the approximate nonlinear equations are modified to account for damping of flywheel rotors, and 
presented in Sec. 2.7 together with all other differential equations employed in numerical simulations 
of nonlinear motion. 

2.1 Derivation of Dynamical Equations 

The system of interest, S, is composed of a rigid body B moving in an inertial or Newtonian 
reference frame N, and several rigid axisymmetric rotors R±, . . . , R p whose mass centers are fixed 
in B. A subset of the rotors R±, . . . , Rp have spin axes fixed in B so that these rotors represent non- 
gimbaled flywheels or reaction wheels. Each of the remaining rotors Rp+ 1 , ... , R p are attached 
to B with one or more massless gimbals which permit the direction of the spin axis to change 
relative to R; these rotors thus represent a number of CMGs, C = p — T . (The latter subset could 
contain gimbaled flywheel rotors as well as CMGs, but we concern ourselves in this work only with 
non-gimbaled flywheels.) This system is illustrated in Fig. 2.1, with rotors R 2 , . . . , R p ~ 1 omitted 
for the sake of clarity. 
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Figure 2.1: Spacecraft with Flywheel and CMG 


The equations of motion are derived using Kane’s equations (Eqs. (6.1.2), Ref. [17]) 

F r + F* = 0 (r = 1, . . . , n) (2.1) 

where F r are generalized active forces for S in N, F * are generalized inertia forces for S in N, and 
n is the number of degrees of freedom of S in N. 

2.1.1 Generalized Speeds 

The system S is holonomic and therefore a complete description of the motion of S in N requires 
n generalized speeds u\, . . . ,u n , conveniently chosen as follows. 

B = uibi + rt 2 b 2 + u 3 b 3 (2.2) 

where N u> B is the angular velocity of B in N, and bi, b 2 , and b 3 are a set of mutually orthogonal, 
right-handed unit vectors fixed in B. 

We introduce unit vectors 0 i fixed in B such that they are each parallel to the spin axis of 
a flywheel rotor Ri, and therefore to the angular velocity B iv Ri of Ri in B. Generalized speeds 
Ui, . . . , ujr + 3 associated with the flywheels can then be used to write the angular velocities as 

B u; R ' = u i+ 3 ft (i = l,...,F) (2.3) 

The inner gimbal of each CMG is fastened to B with a revolute joint in a single-gimbal con- 
figuration, whereas a double-gimbal CMG has the inner gimbal attached with a revolute joint to 
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an outer girnbal, which is in turn mounted in B with a second revolute joint. The axis of each 
revolute joint is assumed to pass through the mass center of the rotor, which is thus fixed in B. In 
order to maintain generality with regard to the number of CMG rotors and gimbals, the angular 
velocities b u) R t+1 , . . . , B u: Rp of the rotors relative to B are not written explicitly, but they must 
be functions of the generalized speeds u^+ 4 , . . . , u n - 3 , where one generalized speed is required for 
every girnbal. 

The final three generalized speeds are associated with the velocity of the mass center S * of S 
in N, 

N v s * = u n - 2 hi + u n -ih 2 + u n n 3 (2.4) 

where hi, n 2 , and 113 are a set of mutually orthogonal, right-handed unit vectors fixed in N. 


2.1.2 Generalized Active Forces 


Let <7 be the set of forces exerted on S except those exerted by B and R 4 on each other (i = 1, . . . ,p). 
The forces in a acting on B,R\, . . . , R p are equivalent to single forces F#, Fi, . . . , F p applied at 
the mass centers B* ,R\, . . . , R* of bodies B,R\, . . . ,R p , respectively, together with couples whose 
torques are Mg, Mi, . . . , M p . 

To account for the forces exerted by B on Ri, we regard them as equivalent to a single force 
F B / Ri applied at R*, together with a couple whose torque is M b / r \ Since R* is fixed in B, F B / Ri 
contributes nothing to F r (i = 1, . . . , p; r = 1, . . . , n). 

Body B makes a contribution ( F r ) B to the rth generalized active force, given by Eq. (4.6.1) 
of Ref. [17]. Likewise, each rotor makes a contribution denoted by ( F r ) Bi . Accounting for con- 
tributions from all bodies in the system, the law of action and reaction, and discounting the 
noncontributing forces, 


F r = (F r ) B + j~](Fr)Rj 

i= 1 

= A r B * • F b + N u R ■ (mb-£ M B/Ri 


1=1 


+ 


1=1 


N Vr R * ■ Fj + N u Ri ■ ( Mi + M B/Ri 


(r = 1 , . . . , n) 


(2.5) 


where N v B * is known as the rth partial velocity of B* in N , N uj b is the rth partial angular 
velocity of B in N, and so forth. 

Since the mass centers of all bodies belonging to S are fixed in B, the system mass center 5* is 
also fixed in B, and the velocities of B * and R* in N can be written in terms of the velocity of S* 
in A, 


jv v b* 

N v R ! 


^v v s* + n u b xr S*B* 
N v s * + N u B x r s * R i 


(* = ,P) 


(2.6) 

(2.7) 


where r s * B * is the position vector from S* to B*, and r s * R i is the position vector from S * to R *. 
Thus, the partial velocities can be written as 


= N v r S * + N u B x 
N v r V = N v r s * + N u B x 
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(r = 1 , . . . , n) 

(* = !,••• r = 1,... ,n) 


(2.8) 

(2.9) 


The angular velocities in N of the rotors can be written in terms of the angular velocity of B 
in N, 


N u R ' = N u B + B u R ' = 


and therefore their partial angular velocities can be written as 


N u Ri = N U3 B + B u Ri (i = 1, . . . , p; r = 1 , . . . , n) 


(2.10) 


( 2 . 11 ) 


Substitution from Eqs. (2.11), (2.9), and (2.8) into (2.5) yields 


F r = ( + N u>? x r s * B * 


F B + N u B ■ l M b + Y Mi 


2=1 


+ [(-vf + N u B x r s *^) . F, ; + B u R ‘ ■ (Mi + M b / r >) 


2—1 


N v? ■ (Fs + ^Fij + 

+ Y BuJ r il • ( M * + M B/Ri 


N u B 


M b + r s * s * xF B + ^ [Mi + r s * R i x F 


2 — 1 


2=1 


= A v r s * • F + N u> B • M + Y2, B “ Ri • [Mi + M B / Ri J (r = 1, . . . , n) 

2=1 

where F is the resultant of the forces in a acting on S, 

p 

F = F B + Y Fi 

2=1 

and M is the moment of a about S *, given by 

p , 

M = M s + r s * B * xF B + ^ (M* + r s *^ x F* 


(2.12) 


(2.13) 


2=1 


(2.14) 


2.1.3 Generalized Inertia Forces 


All bodies belonging to S contribute to the rth generalized inertia force F* according to Eqs. (4.11.5)- 
(4.11.7) of Ref. [17]. As pointed out in part (c) of Problem 8.18, Ref. [17], the inertia torque in N of 
a rigid body B is simply the negative of the time derivative in N of the central angular momentum 
of B in N, hence, 

f; = (f;)b + Y( f ;) Ri 


2=1 


N 


d 


_ I N B* , rrl N a B*,N,,B u NttB/B* 
— — I v r • m B a + cJ r • — ' 

P / N , 


Y [ N v r Rt ■ rrii N a R i + N u r R > • | W ) (r = 1, . . . ,n) 


2=1 


(2.15) 


where tub is the mass of B, N a B * is the acceleration of B * in N, and so forth. Substitution from 
Eqs. (2.11), (2.9), and (2.8) into (2.15) yields, after some manipulation, 

/ p \ p N 


F* = - J V 


N„S* 


m B N a B * + ^ 
^ 2=1 


m; N a - 


V B u Ri ■ — N n Ri/R 

r rtf 


i = 1 


- N u B 


'N 


dt 

N 


« N U B/B* + r S*B+ xmB N a B* + Y « N u R i/R ! + f 
nt ' \ nt 


i= 1 


s * R i x rrii N a R t 


(r = 1, . . . , n) (2.16) 


After noting once more that 5 1 *, H*, and i?*, . . . , R* are all fixed in B, the accelerations of B* 


and R* in N can be written in terms of the acceleration of S* in N, 

N a B * = N a s * + N a B x r s * B * + N u B x N u> B xr s * B * 
N a R i = N a s * + N ot B x x V x 

where is the angular acceleration of B in N. Thus, 

r s * B * x m B N a B * + £ r s * R * x m, ^ 

2=1 

= r s * s * x m B ( N a s * + N a B x r s * B * + N u B x N u B x r s * B *) 


(2.17) 

(i = l,..., / 9) (2.18) 


2=1 


+ x mj( ^ a s * + x + N u, B x N u> B x r s * R *) 

+ £m i r s **A x N a s * 

9 1 

+ £ 


= m^r 




+ <m B 


2=1 




U - r S*B* r S*B* 


mi 


2=1 




+ N u> B x < m B 


= 0 x N a s * 


V i= i 


,S*B * 


R*/S* 


U _ t s*b* t s*b* 


+ £ 


mi 


i = 1 


U -r 




S*R* r S*R* 


N a B 


U - r s * R t r s * R i 


N a B + N U B xi I B */ S * + . iV w . 




(2.19) 


2=1 


where the vector that premultiplies 


TV „ S* 


vanishes because of the definition of the system mass 


center, and where U represents the unit dyadic. As described in Sec. 3.6 of Ref. [17], I B ^ denotes 
the inertia dyadic relative to S* of a fictitious particle situated at B* and having a mass equal to 
m B . A similar description applies to the dyadic I R ^ s . 

The central angular momentum of B in N is simply 


N 


U B /B* = jB/B* . 


B 


(2.20) 


where 1 B / B * is the inertia dyadic of B relative to B*. Hence the time derivative of B / B * in N 

is given by 


N 


d_ JV R b/b * = jB/B* .N a B + N u B x j B/B* . N „ B 

dt 


(2.21) 
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The central angular momentum of any rotor in N can be written as 

N-^Ri/R* _ ^Ri/R* ' ( N u B _|_ B u R i ) 

= I Ri/R i • N u B + B n Ri/R i (i = l,...,p) (2.22) 

where l Ri / R * is the inertia dyadic of Ri relative to R*, and B H Ri ^ is called the central angular 
momentum of in B. The time derivative of N H Ri / R i in the case of a flywheel rotor is different 
from that of a CMG because the mass distribution of the rotor does not change relative to B in 
the former case, whereas it does change in the latter. For flywheel rotors, 


N 


^N H Ri/R* 

dt 


B 


= ^LN-ftRi/R* + N^B x N H Ri/R* 

dt 


B 


= 1 fl Ri/ R t . N U B + B H Ri/Ri} + N UJ B x (j Ri/R* . N^B + B^Ri/R* 


B 


d 


= 0 • N u B + 1 R i/ R i . N ot B + N u B X I Ri / R i ■ N U B + — B H R i/ R i + X B HW 

dt 

(i = l,...,R) (2.23) 


where the zero dyadic 0 results because there is no change in the inertia dyadic l Ri / R i with respect 
to time in reference frame B. Unlike the flywheel rotor, the CMG rotor’s inertia dyadic does change 
with respect to time in B, and in view of Problem 5.14, Ref. [17], 


Bj Ri, 

J Ri/R* — _J Ri/R* I B u> R i X I R i/ R i _ J R i/ R i x B (_J R i 

dt~ dt~ 

= 0+ B LL, Ri xI Ri/R i -l Ri/R i x B uj Ri (i = R+l,... ,p) (2.24) 


This result is then used in place of 0 in Eqs. (2.23) and yields 


N 


d_ 

dt 


B 


— N-^Ri/R* = j Ri/R* . N^B + N^B x j Ri/R* . N^B + ± B-^Ri/R* + N^B x fljjfli/flf 

dt 


+ ( B u Ri x I Ri / R i -l R i/ R t x B u> Ri ^J ■ N u B (* = F+l,... ,p) (2.25) 


After noting that mg N a B * + X)( , =1 ra; a R t = mg Af a s *, where ms is the mass of system S 
and N a‘ s * is the acceleration of S* in N, and then substituting from Eqs. (2.25), (2.23), (2.21), 
and (2.19) into (2.16), one obtains 


F* = 


.m s N a s *- N u r B 


P / B 


I S/S* . 


N„,B , N, . B 


' a “ + iy u ” x I 5 / 5 ■ “ U 


N , , B 


+ J2 ( j t B n Ri / R i + N u B x B H Ri / R * 


i = 1 
P 


+ ^2 {^ B UJ Ri X l Ri / R * —J. Ri / R i x B CJ Ri ^j • N U) B 
i=F + 1 

P n j 

-E 


B, .Ri 


i=l 


'u;" 1 - — N H Ri / R i (r = 1, . . . , n) 


(2.26) 
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where I 5 ’/ 15 ' , the inertia dyadic of S for S*, is given by 


I S/S* = jB/B * + J B*/S* + f I Ri/Rt + I R-/S 


1=1 


(2.27) 


2.1.4 General Vector-Dyadic Equations 

According to Eqs. (2.1) the generalized inertia forces from Eqs. (2.26) can now be added to the 
generalized active forces from Eqs. (2.12) to yield vector-dyadic equations of motion for a spacecraft 
containing flywheels and CMGs 


7V v S* 


+ N v B 


F — m s N a 5 * 

M - 


'B 


d 


I S/s* . N a B + N^B xl S/S* . N^B 

+ J2[j t B K IU/R !+ N u> ^ 

1=1 \ 

+ ( B UJ Ri xI Ri ^ R i — I Ri / R i x B U) Ri j • N u 


■ B X Ri / R i 


i=T + 1 
p 


. B 


N 


+£ 

2 — 1 


B U) Ri 


d 


Mi + M B/Ri - — N H R i/ R i 1=0 (r = 1, . . . , n) 


(2.28) 


Eqs. (2.28) are completely general with regard to the number and orientation of flywheel rotors, 
and the number of CMG rotors and gimbals. In this form, they are applicable to variable-speed 
CMG rotors. These equations of motion, and the expression for generalized inertia forces, can be 
compared to previous work. 

Reference [18] is concerned with gyrostats and relevant equations that can be dealt with easily 
by an analyst, and quickly by a computer. Expressions for generalized inertia forces are presented 
separately for a gyrostat containing a single cylindrical rotor, and for one containing a single 
spherical rotor; an underlying general relationship (C65) developed in Appendix C can be shown 
to give rise to Eqs. (2.26) presented here, when no CMGs are present {p = J 7 ). 

The term in the second line of Eqs. (C65) accounts for but a single rotor, although additional 
rotors can be handled straightforwardly by adding a sum of similar terms. By appealing to Eqs. (24), 
(C61), and (C35) of Ref. [18], and replacing their labels G, B , and A for the gyrostat, rotor, and 
carrier respectively with our S, Ri , and B, we can write 

p 


—F* = N v r s * . N F S + N u B NrrS 


t s +E 


B u R ' 


i = 1 


Ri/Ri . N Ot Ri + N UJ Ri X I Ri/ R i . ^U} Ri ^j 


N . 


N 


NS* 

= v r • ms a 


d 


N » s * ■ N u B • ^H S / S * + ^V. “"nRi/V (r = 1 , . . . , n) 

( JjL lll 


2—1 


(2.29) 

The central angular momentum in N of a gyrostat containing a single rotor, expressed by Eqs. (C24) 
and (3) of Ref. [18], is easily modified to include multiple rotors as in Eq. (14) in Sec. 6.1 of Ref. [19], 

W H S / S * = I S / S * • N u B + J2l Ri / R i • B U Ri = I S / S * .Vlf B^Ri/R* ( 2-30 ) 


2—1 


1=1 
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and differentiation with respect to time in N yields, after reference to Eq. (C34) of Ref. [18], 


v J 

/± S/S * _ J S/S* m N B 

dt 


N 


cj b xI s / s * . N u> B + J2 ( ^- t B H Ri/R * + n lj b x B H Ri / R i 


i=l 


(2.31) 


Substitution from Eq. (2.31) into (2.29) clearly reproduces Eqs. (2.26) when CMGs are absent. 

The use of the system mass and inertia scalars, together with the moment of inertia for the 
axis of symmetry of each rotor, is shown in Ref. [18] to lead to greater efficiency than use of mass 
properties of individual bodies in a gyrostat; this advantage happens to accrue to Eqs. (2.26) and 
(2.28) developed here. Additionally, in the interest of efficient equations of motion, the authors of 
Ref. [18] advocate the use of the partial angular velocities N u Ri of a gyrostat rotor in N , rather 
than partial angular velocities B u Ri in B, when motion of the rotor relative to B is not prescribed 
(known a priori). The latter are chosen here because they bring into evidence S H Ri / R i and its 
time derivative in B, quantities which are needed for the forthcoming comparison in Sec. 2.3 to 
relationships used by others for dealing with CMGs. 

Rheinfurth and Carroll present in Ref. [20] vector-dyadic equations of motion (27) for a space- 
craft composed of a rigid carrier and a rigid appendage whose mass center is fixed in the carrier. 
Additional appendages are accounted for easily by forming a sum, as they do in Eq. (16), but the 
resulting vector-dyadic expression will give rise to only three scalar equations. It is pointed out 
near the bottom of p. 6 that a CMG can be regarded as an appendage; however, the motion of 
every appendage relative to the carrier must be prescribed if the three relationships are to serve as 
dynamical equations governing the motion of the carrier. Reference [20] does not contain counter- 
parts to the n — 6 of our Eqs. (2.28) that govern the motion of the rotors, or to the three equations 
governing the translational motion of the system. It is now shown that Rheinfurth and Carroll’s 
Eqs. (27) give way to the first three of Eqs. (2.28) here when all rotors are permitted to be CMGs 
(J- = 0). After forming the required sum, and replacing their symbols L with M, I with I 5 / 5 , I p 
with , O with N u: B , lu p with B u> Ri , (0) v with N a B , and (u; p ) p with B a Ri , one obtains 


M = I s / s * • N a B + N u B x • N u B + 


i= 1 l 


J Ri/R* , B qRi _|_ B -Ri ^ j Ri/R* . B^Ri 


+ ^ B U) Ri X I Ri / R i — X R i! R i X B U) Ri ^ • N (jl) B + N U) B X \ Ri / R i • B U Ri 


(2.32) 


After noting that l Ri / R i ■ b uj R >- = and I Ri / R i • B a R% + B u> Ri x l Ri / R i . B u> Rl = 

B d B H R ^ R i /dt, Eq. (2.32) can be rewritten as 

M = I s / 5 * • N a B + N u B x l s / s * • N u B (2.33) 


+ E 


i— 1 _ 


B 


/L B fX TR / R * + N (ji) B X B H Ri / R i + ( B U) Ri X I Ri / R i - I R i/ R i x B u> Ri ) • N UJ B 


It is evident that forming dot products with this expression and three vectors n cj b produces the 
first three scalar relationships given by Eqs. (2.28). 


2.2 Dynamical Equations for A Complex Gyrostat 

A spacecraft known as a simple gyrostat is described in Sec. 3.6 of Ref. [16]; the system S in the 
preceding discussion becomes a simple gyrostat when the number of flywheels /F is equal to 1, and 
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Figure 2.2: Spacecraft with Flywheels 


when no CMGs are present (C = 0, thus p = J-). A spacecraft with more than one flywheel, such 
as the one shown in Fig. 2.2, will be referred to as a complex gyrostat; equations of motion with 
T = 6 are given in the following material, and are shown to correspond to those given in Sec. 3.7 
of Ref. [16] when T = 1. 

Without a great loss of generality one can at this point work with six flywheel rotors R\, . . . ,Rq 
(R = 6 ) arranged in three counter-rotating pairs as shown in Fig. 2.2, with the spin axes of R\ and 
7?4 parallel to bi, Rg and i ?5 parallel to b» 2 , and R 3 and Rq parallel to b> 3 . The generalized speeds 
114 , ... , ug associated with the flywheels are then used to form the angular velocities B u Ri of Ri 
in B, (i = 1, . . . , 6 ) 


Rl = U4 bi, 

B UJ R2 = u 6 b 2 , 

B u> Rs = u 8 b 3 , 

' R4 = U 5 bi, 

B w fis = U 7 b2, 

B U) R& = Ug b 3 


(2.34) 


Since the axis of symmetry of each flywheel remains fixed in B , and is in the case of R\ parallel 
to bi, the central inertia dyadic of R\ can be expressed as 

iRi/Rl = jt 1 bi + 7r(b 2 b 2 + b 3 b 3 ) = Jbibi + K { U - 6464 ) (2.35) 


where J and K are the central principal moments of inertia of a flywheel rotor for, respectively, the 
axis of symmetry, and any line perpendicular to the axis of symmetry. In fact, the inertia dyadics 
for the 6 flywheel rotors are given by 

l Ri / R * = l R '+3/ R l+3 = j^b* + K { U - bibi) (i = 1, 2, 3) (2.36) 
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The angular velocities of the rotors in B are 


B ui Ri = 


U2i+2^>i 

U2i-3^>i-3 


(*=1,2,3) 
(i = 4,5,6) 


(2.37) 


In order to evaluate Eqs. (2.28) one must develop expressions for angular momenta of Ri ap- 
pearing there, and time derivatives and sums thereof. In view of Eqs. (2.34), (2.2), and (2.10), 


Bj^Ri/R* _ j Ri/R* m B^Ri _ 


Jbjbi + If (U — bjbj) • U2i+2^>i (* = 1,2,3) 

■7t»j_ 3 bj_ 3 + If (U — bj_ 3 bj_ 3 ) • u 2 j- 3 b ?; _ 3 (* = 4,5,6) 


Ju2i+2^i 

JU2i-3^>i-3 


(* = 1,2,3) 
(* = 4,5,6) 


(2.38) 


E B H Ri/R * i = J[(**4 + U 5 )bi + (*i 6 + U 7 )t>2 + (U 8 + Ug)b 3 ] 


i = 1 


(2.39) 


N u> B x 


E* H 


Ri/R* = 


i = 1 


(2.40) 


J{[u 2 (u 8 + u 9 ) - u 3 (u 6 + u 7 )]bi + [U3(U4 + U 5 ) - ui(u 8 + **g)]b 2 + [**i(?*6 + u 7 ) - ** 2 (**4 + u 5 )]b 3 } 


B ** 

A B-ftRi/R* _ f J^2i+2^i (* = 1,2,3) 

dt \ J tt2i— 3 bj_ 3 (* = 4,5,6) 


(2.41) 


6 Bj 

E ^ = J[(**4 + *t 5 )bi + (u 6 + ** 7 ) 1*2 + (**8 + **9)b 3 ] (2.42) 


N-^Ri/R* 


j Ri/R* m N ^ Ri 


Jb 4 bi + If (U - bibi) 

• NfjjRi 

>Ibj_ 3 bj_ 3 + K (U - b 

i-3bj- 3 ) 


f (J — K){ui + tt 2 i+ 2 )bj + K N u Ri 

\ (J - K)(Ui-3 + **2i-3)t*i- 3 + K N U} Ri 


(1 = 1,2,3) 

(* = 4,5,6) 

(* = 1,2,3) 

(* = 4,5,6) 


(2.43) 


N, Ri, 

_ A r H Ri/R* = ±Nj^Ri/R* + N u Ri x *V Ri/R* 
dt dt 


= '£ (l^/ R i . N u R ^ + N u Ri x (l Ri / R * ■ N u R ^ 

= 0 • N uj Ri + I^ ,; /*** . ^ a R i _|_ N ^Ri x . N^Ri 


(2.44) 
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but 


N, B , 

N a R, = ±N u R t= * N Ri + N u}Bx N u R i 

dt dt 


B 


d 


= j t ( N u B + B u R *)+ N u B x( N u B + B u Ri ) 
B 


N a B 


± B u R i+ N u> B x B u R " 
dt 


N a B + «2i+2bj + N w B x U2i+2bj 
N a B + «2i_3bi_3 + N U> B X U2iSi-2 


(* = 1,2,3) 
{i = 4,5,6) 


(2.45) 


so 


jRi/R* . N a R i = 


Jbibi + if (U - bibi) 

• N a Ri 

Jbi-sbi-g + if (U - b 

i-3bi- 3 ) 


N Ot Ri 


(* = 1,2,3) 

(* = 4,5,6) 

jf Jv a B -b i + U2i +2 + 0j bj + if(U - bib*) • N a Ri 
J ( N a B • bj_ 3 + u 2 i-3 + O) bj_ 3 + if (U — bj_ 3 bj_ 3 J ■ 'a 
J (iti + U2i+2) hi + K (U - hihi) • N a Ri (i = 1,2,3) 


N^Ri 


and 

N UJ Ri X \ Ri / R *i . N U> Ri = 


J (iii - 3 + ii 2i _ 3 ) bj-3 + K ( U - bj_ 3 bj_ 3 ) • N a Ri (i = 4,5,6) 

(* = 1,2,3) 


(2.46) 


N u> Ri x 


N U) Ri X 


(J-K)( Ui + n 2 i+ 2 )bi + if ^ 

(J - K)(v,i- 3 + ** 2 i— 3 )bj — 3 + K N u Ri (i = 4,5,6) 

(J - K)(ui + « 2 i+ 2 ) xbj (i= 1,2,3) 

(J - if)(ui _ 3 + U 2 i- 3 ) N u Ri x bj _3 (i = 4,5,6) 


(2.47) 


Therefore, substituting from Eqs. (2.47) and (2.46) into (2.44), one obtains 

N, 

_ Nj^Ri/R* = 

dt 


(2.48) 


J(ui + ii 2 i+ 2 ) b* + if(U - bjbi) • N cx Ri + ( J - K)(m + U 2 i +2 ) N u Ri x b* 

(* = 1,2,3) ^ 

J (ui -3 + ii 2 i— 3 ) bj -3 + if (U - bi_ 3 bj_ 3 ) • N a Ri + ( J — K)(ui-s + 1121 - 3 ) x bj_ 3 
(* = 4,5,6) 

Now write M as 


M = M\b\ + M 2 b 2 + M 3 b»3 


(2.49) 


and let unit vectors bi, b 2 , and b 3 be parallel to central principal axes of inertia of S, so that 

jS/s _ + i 2 b 2 b 2 + 13 ( 331,3 (2.50) 

where I\, i 2 , and I 3 are central principal moments of inertia of S. 

By inspecting Eqs. (2.2), (2.34), and (2.4), one can construct a table of partial angular velocities 
and partial velocities necessary for forming scalar relationships from Eqs. (2.28). Upon referring to 
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Table 2.1: Partial Angular Velocities and Partial Velocities for Complex Gyrostat 


r 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

N u B 

bi 

b 2 

b 3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

B u Rl 

0 

0 

0 

bi 

0 

0 

0 

0 

0 

0 

0 

0 

B U J R2 

0 

0 

0 

0 

0 

b 2 

0 

0 

0 

0 

0 

0 

B <jl> Rs 

0 

0 

0 

0 

0 

0 

0 

b 3 

0 

0 

0 

0 

B U) Ri 

0 

0 

0 

0 

bi 

0 

0 

0 

0 

0 

0 

0 

B U) R 5 

0 

0 

0 

0 

0 

0 

b 2 

0 

0 

0 

0 

0 

B u Re 

0 

0 

0 

0 

0 

0 

0 

0 

b 3 

0 

0 

0 

N S* 

V r 

0 

0 

0 

0 

0 

0 

0 

0 

0 

ni 

n 2 

n 3 


Table 2.1, it is evident that the terms involving N v B * and B u: Ri (i = 1, . . . ,6) contribute nothing 


to equations of motion (2.28) for r = 1, 2, and 3. Substituting from Eqs. (2.50), (2.49), (2.42), and 
(2.40) into (2.28), one obtains 

hui + J(ii 4 + U 5 ) = {h - h)u 2 Uz - J[u 2 {us + u 9 ) - u 3 (u 6 + u 7 )] + Mi (2.51) 

I 2 U 2 + J{iiQ + ii 7 ) = (h ~ h)uiu 3 - J[u 3 (u A + u 5 ) - ui(u 8 + w 9 )] + M 2 (2.52) 

I 3 U 3 + J(ii 8 + U 9 ) = (h - h)uiu 2 - J[ui(u 6 + u 7 ) - u 2 (u4 + u 5 )] + M 3 (2.53) 

It is also apparent from Table 2.1 that N v B * and N u: B contribute nothing to Eqs. (2.28) for 

r = 4, . . . ,9, and substitution from Eqs. (2.48) into (2.28) yields 

J{iii + u A ) = (Mi + M B/Rl ^j ■ bi (2.54) 

J(iii + u 5 ) = (m 4 + M s/i?4 ) • bi (2.55) 

J{ii2 + u 6 ) = (M 2 + M b ^ • b 2 (2.56) 

J(ii 2 + 117 ) = (m 5 + • b 2 (2.57) 

J(ii 3 + u 8 ) = (m 3 + M b ' r 3 ) • b 3 (2.58) 

J{ii 3 + ii 9 ) = (M 6 + M b / r 6 ) • b 3 (2.59) 


Although translational equations of motion are not used in what follows, they are obtained easily 
from differentiation of Eqs. (2.4) with respect to time in N, and from Eqs. (2.28) with r = 10, 11, 
and 12, and are simply 


msu r = F r ~ 9 (r = 10, 11, 12) 

where 


F = Eihi + F 2 n 2 + F 3 n 3 


(2.60) 


(2.61) 
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Eqs. (2.51)-(2.60) are thus a complete set of nonlinear dynamical equations of motion for a complex 
gyrostat composed of a body B and six axisymmetric rotors Ri , . . . , Rg whose mass centers and 
spin axes are fixed in B\ the rotors are arranged in pairs, with the spin axes in each pair parallel 
to each other and to a central principal axis of inertia of the gyrostat. 

Eqs. (2.51)-(2.59) exhibit dynamic coupling; numerical solution of these differential equations 
can proceed more efficiently when the equations are decoupled, and this is accomplished easily by 
substituting from, for example, Eqs. (2.54) and (2.55) into (2.51). Consequently, Eqs. (2.51)-(2.53) 
give way to 

(h ~ 2J)iii = (I 2 - I 3 )u 2 u 3 - J[u 2 (u 8 + u 9 ) - u 3 (u 6 + u 7 )] + Ml 

- (Mi + M B/Rl + M 4 + M b/R4 ^J • bi (2.62) 

(h - 2 J)u 2 = (h - h)u\u 3 - J[u 3 (ui + u 5 ) - ui(w 8 + u 9 )] + M 2 

- (m 2 + M b/R2 + M 5 + M B/R ^ • b 2 (2.63) 

(I 3 - 2 J)u 3 = ( h - I 2 )uiU 2 - J[ui(ug + u 7 ) - 112(114 + w 5 )] + M 3 

- (m 3 + M B / Ri + M 6 + • b 3 (2.64) 

which yield hi, h 2 , and u 3 effortlessly and, with these in hand, h 4 , . . . , iig are obtained straightfor- 
wardly from Eqs. (2.54)-(2.59). 

In Sec. 3.7 of Ref. [16], Kane, Likins, and Levinson give rotational equations of motion for a 
simple gyrostat, where their body A plays the part of B in the preceding discussion. Eqs. (3.7.28)- 
(3.7.31) are obtained by supposing the spin axis of the rotor to be parallel to the third central 
principal axis of inertia of the gyrostat, therefore Kane’s unit vector (3 plays the part of b 3 , and 
B plays the part of rotor R 3 (or Rq). To show the correspondence between the aforementioned 
equations and Eqs. (2.51)— (2.59), note that u>i, cu 2 , and w 3 play the parts of u 4 , u 2 , and u 3 respec- 
tively, that a oj b plays the part of u 8 , and set generalized speeds u 4 , . . . ,u 7 and ug and their time 
derivatives to zero. Eqs. (2.51), (2.52), and (2.53) then give way to 


hiii = (I 2 - I 3 )u 2 u 3 - Ju 2 u 8 + Mi 

(2.65) 

I 2 U 2 = (I 3 ~ h)uiu 3 + Ju\u 8 + M 2 

(2.66) 

I 3 u 3 = (h — h)uiu 2 — Ju 8 + M 3 

(2.67) 


which are identical to Eqs. (3.7.29)-(3.7.31) of Ref. [16]. In the absence of rotors R\, R 2 , R 4 , R 5 , 
and Rq, Eq. (2.58) is the only one of Eqs. (2.54)-(2.59) that remains applicable. Solving Eq. (2.67) 
for ii 3 and then substituting from the result into (2.58), one obtains 

1 - i.) jjg = ( M3 + MJi/Ril ) ’ _ Vi - h)mu 2 + M 3 

h) J h 

which is the same as Eq. (3.7.28) in Ref. [16] provided the external torque M 3 acting on the rotor 
vanishes, or, to be more precise, the dot product M 3 • b 3 vanishes. The disappearance of this 
dot product is perfectly reasonable if the rotor does not possess a magnetic dipole moment and 
is housed inside B where it is protected from the action of aerodynamic forces, in which case the 
principal contribution to M 3 is gravitational moment exerted by the celestial body about which 
the gyrostat orbits; an expression for an often used approximation of gravitational moment is given 
in Eq. (2.6.8) of Ref. [16], where it becomes clear that the dot product in question vanishes because 
b 3 is parallel to an axis of symmetry of the rotor. 
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Figure 2.3: Spacecraft with CMG 


2.3 Dynamical Equations for Spacecraft with CMGs 

In Sec. 2.2, Eqs. (2.28) are tailored to describe the motion of a spacecraft carrying flywheels and no 
CMGs; in this section, the important case of a spacecraft with CMGs and no flywheels is examined 
(J 7 = 0, p = C). In Ref. [6], Wie et al. develop a scheme for controlling a spacecraft’s attitude 
and managing the angular momentum of a collection of CMGs, and the results are applied to a 
space station. Under certain conditions Eqs. (2.28) can be shown to give rise to the relationships 
employed in Ref. [6] as a basis for control law design. 

One may recall that Eqs. (2.23) are developed differently from Eqs. (2.25) because the central 
inertia dyadic of a flywheel rotor does not change with respect to time in B whereas that of a CMG 
rotor does; consequently, the second sum in Eqs. (2.28) is made up entirely of contributions from 
CMGs but not from flywheels. Each CMG rotor Ri is attached to an inner gimbal Gi (treated as 
massless) as shown in Fig. 2.3; the axis of symmetry of Ri remains fixed in Gi, and the relative 
angular speed G of Ri in Gi remains constant. The angular velocity of Ri in B can thus be expressed 
as 


B u Ri = B u Gi + Gi u Ri = B u Gi +n*/i (i = JF+l,... ,p) (2.69) 

where qq is a unit vector parallel to the axis of symmetry of Ri, and the angular velocity B u: Gi of 
the inner gimbal in B accounts for the motion relative to Gi and to B of an outer gimbal (shown, 
but not labeled in Fig. 2.3) if one is present. The central inertia dyadics of the CMG rotors may 
be expressed as 


l Ri/R * = Jlili + £(U - 7j7j) (i = T + 1, • ■ ■ , P) (2.70) 
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where J and /C are the central principal moments of inertia of a CMG rotor for, respectively, the 
axis of symmetry, and any line perpendicular to the axis of symmetry. With the aid of Eqs. (2.69) 
and (2.70), the second sum in Eqs. (2.28) can be written as 

p 

'y ^ U} Ri ^ \ R i/ R *i — \ R i/ R t X B (jj 

i=T + 1 

= y B “ Ri } 

i=T + 1 

= E u-nxV)] 

i=T + 1 

= E [{ B u Gi + Qji) x yfii -7i7i x ( B u Gi +O7J] +0} 

i=T + 1 

= E Xlili-lili X B ^ Gi ) (2.71) 

i=J 7 +l 


where the zero dyadic 0 results because, for any vector v, it can be shown that v x U = U x v. 

Now, this second sum can be neglected in comparison to the first sum after considering typical 
numerical values involved in each. According to Ref. [21], p. 2.3-5, the International Space Station 
CMG rotor spin speed 11 is constant, 6,600 rev/min or 691 rad/s. On p. 2.3-7, typical girnbal speeds 
are reported to be 3.5xl0 ~ 3 rad/s (0.2 deg/s), with maximum girnbal speeds of 5.4xl0 -2 rad/s 
(3.1 deg/s), which gives an idea of the magnitude of B u Gi . The central angular momentum of Ri 
in Gj, g *H Ri / R * = I R i/ R * . Gi u ) Ri , has a magnitude of JSl = 4,745 N-rn-s (3,500 ft-lbf-s, often 
referred to as “the” momentum of a CMG) which is approximately the magnitude of S H Ri / R i 
since the girnbal speeds are so much less than the spin speed. From the values of momentum and 
spin speed, one may deduce that J = 6.87 kg-m 2 . When one compares the magnitude of 
that appears in the first sum to the product of J and an approximate magnitude of B m Gi that 
appears in the second sum, one finds that the former, 4,745 N-m-s, is some 4 orders of magnitude 
larger than the latter, (6.87 kg-m 2 )(5.4xl0 ~ 2 rad/s) ~ 0.4 N-m-s, and neglecting the second sum 
is easily justified. 

The first sum in Eqs. (2.28) can be denoted by the vector r, 


SB 


T = E ( j t B n Ri/R * + N u B X B n Ri / R t 


i=F + 1 
B 


d_ 

dt 


Y b hW 

^=JE-+l J 


B 


Y^ B H R i/ R i 


(2.72) 


where the term in parentheses, the resultant of the central angular momenta in B of all CMGs, 
can be defined as the vector h 

h= Y B tl Ri/R i (2.73) 

i=B + 1 

and the time derivative of h in B is thus 

B ■, 

h = - 7 -h (2.74) 

dt v ' 


19 


Substitution from Eqs. (2.74) and (2.73) into (2.72) yields 

r = h + N u B x h 


(2.75) 


which can be recognized as a vector form of Eqs. (3) in Ref. [6], three scalar equations that result 
from forming dot products between Eq. (2.75) and unit vectors bi, b2, and b3, fixed in B but not 
necessarily parallel to central principal axes of inertia of S. The measure numbers appearing in 
Eqs. (3) of Ref. [6] are related to the vectors in Eq. (2.75) as follows, h r = h • b r , h r = h • b r , 

u r = N u } B • b r , and u r = r • b r [not to be confused with u r defined in Eq. (2.2) here], (r = 1,2,3). 

Although r represents in some sense the effect of CMGs on B, one must remember that it fails 
to account for those effects in the neglected second sum; furthermore r is not the resultant moment 
of forces exerted by the CMGs on B, for this resultant is given by — by virtue of 

the law of action and reaction. 

In addition to Eqs. (3), the control law design in Ref. [6] is based upon the three scalar Eqs. (1) 
appearing therein, and one may obtain these from Eqs. (2.28) by reasoning in the following manner. 
First, the term involving N v B * contributes nothing for r = 1,2,3, nor does the term containing 
B u> Bi . Although generalized speeds corresponding to angular gimbal speeds have not been intro- 
duced, the partial angular velocities B u) Bi of the CMG rotors in B vanish for r = 1,2,3, just as 
they do in the case of flywheel rotors (see Table 2.1). Next, the second sum is neglected as before 
and substitution is made from Eq. (2.72) into (2.28). Since N u> B = b r , the first three of Eqs. (2.28) 
become 

(l s / s * . N a B + N u B x I 5 / 5 * . N u B + r) • b r = M . b r (r = 1,2,3) (2.76) 

The purpose of control laws such as the one designed in Ref. [6] is to determine what r should be 
in order to control the attitude of B, and minimize the magnitude of h. One may regard M as the 
sum of two terms, the gravitational moment exerted on S, and all other contributions to M. Using 
an approximate expression for gravitational moment such as that given by Eq. (2.6.3) in Ref. [16], 
and the vector w to represent the remaining contributions, Eqs. (2.76) can be rewritten as 

(l s / s *. N a B ) •b r =(-Vxf s *.V + 3r i 2 rxf s *.f-T + w) • b r (r = 1,2,3) 

(2.77) 

which are seen to be identical to Eqs. (1) of Ref. [6] once the measure numbers appearing therein 

are defined as uii = N cx B • b r , u r = n <jJ B • b r , c r = r • b r , u r = r • b r , and w r = w • b r 
(r = 1,2,3). Here the symbol n represents angular speed of a circular orbit, rather than the 
number of degrees of freedom of a system; this will be the meaning of n henceforth when used in 
an equation, unless it is used as the limit of an index such as r, or in the index of a generalized 
speed u or its time derivative. The moments and products of inertia appearing in Eqs. (1) are 

defined as I rs = b r • I 5 / 5 • b s (r, s = 1,2,3); their absolute values range from 0.15x10® to 80xl0 6 
kg-nr 2 in Ref. [6], and are currently expected to be in the range of 1 x 10® to 200x10® kg-rn 2 for the 
Assembly Complete configuration of the International Space Station. Reorientation of a few CMG 
rotors with central inertia scalars on the order of J = 6.87 kg-rn 2 changes the mass distribution of 
S' by a negligible amount; therefore, I rs may be regarded as constants although this is in fact an 
approximation . 

A note is now in order regarding the number of differential equations (2.28) on the one hand, and, 
on the other hand, the number of differential equations employed in the approach of Ref. [6]. The 
number of dynamical equations of motion (2.28) is equal to the number n of degrees of freedom 
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of S in N, exactly the number required to specify the motion of S completely. Three of these 
equations describe the translation of S*, and the remaining n — 3 govern rotational motions of B 
and the CMG rotors. In contrast, rotational motion is represented in Ref. [6] by just 6 differential 
equations — 3 of Eqs. (2.77) and 3 of Eqs. (2.75) — no matter what is the actual value of n. 


2.4 Approximate Equations for a Spacecraft with 
Flywheels and CMGs 


In Sec. 2.2, Eqs. (2.28) are applied to a spacecraft carrying only flywheels, whereas they are applied 
in Sec. 2.3 to a spacecraft containing only CMGs. The results from both Sections are combined 
here in order to produce equations of motion applicable to a spacecraft with both kinds of attitude 
control devices; the equations thus formed must be considered approximate as they include the 
considerations set forth in Sec. 2.3. 

Even though the second sum in Eqs. (2.28) receives no contribution from the flywheel rotors, 
one can deal with their contribution to the first sum by forming counterparts to Eqs. (2.72)— (2.75), 


_ A 
T = 


B 


dt 


T 


Y b \ i R i/ R i 


v i=l 


+ Vx 


T 


Y j B a Ri / R i 


\i=l 


(2.78) 


H = ^ s H Ri/i ^ (2.79) 

i = 1 



H = 

dt 

(2.80) 

and 

f = H+ N u B xH 

(2.81) 

In view of Eqs. (2.39), (2.42), and (2.40), the following definitions can be made: 


Hi = H • bi = J(u 4 + u 5 ), 

Ho = H • b 2 = J{uq + u 7 ), H 3 = H-h 3 = J(u 8 + u 9 ) 

(2.82) 

Hi = H ■ bi = J(w 4 + ft 5 ), 

Ho = H • b 2 = J(v,6 + u 7 ), H 3 = ll-b 3 = J(u 8 + u g ) 

(2.83) 


fi = t • bi = Hi + u 2 H 3 - u 3 H 2 

(2.84) 


t 2 = t • b 2 = H 2 + u 3 Hi - u\H 3 

(2.85) 


t 3 = f • b 3 = H 3 + UlHo, - u 2 Hi 

(2.86) 

Using the same reasoning that preceded Eqs. (2.77), and substituting from Eq. (2.78) as well 
as (2.72), the first three of Eqs. (2.28) become 

N a B ) -b r 
= (-Vxl 5 / 5 *. 

n lo b + 3n 2 r x I s / s • r — r — r + wj • b r (r = 1, 2, 3) 

(2.87) 
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In the term representing the gravitational moment, the coefficient n 2 results from an assumption 
that S* travels in a circular orbit about the mass center P * of a primary celestial body. Unit vector r 
has the direction of the position vector from P * to S*, and is fixed in a local-horizontal- local- vertical 
reference frame L whose angular velocity N u L in N has a magnitude of n. 

In addition to the dynamical equations of motion dealt with thus far, kinematical differential 
equations for variables describing the orientation of B in some reference frame are required; it is 
convenient to choose L as this frame and to employ angles belonging to a body-three, 2-3-1 rotation 
sequence as set forth on p. 423 of Ref. [16]. The sequence is also known as pitch- yaw-roll, with the 
angles denoted by 8 1, 0 2 , and 83 respectively. An associated direction cosine matrix is given by 



li 

h 

h 


bi b 2 

CiC 2 S 3 S 1 - C1S2C3 
s 2 C 2 C 3 

—S\C 2 SiSaCs + SsCi 


t»3 

C\S 2 S 3 + C 3 S\ 

-c 2 s 3 

C3C1 - SiS 2 S 3 


( 2 . 88 ) 


where Si = sin(9i, C 2 = cos0 2 , and so forth, and where li, 1 2 , and I3 are a right-handed, mutually 
perpendicular set of unit vectors fixed in L, with I3 = — r. n uj l is parallel to 1 2 , 


n uj l = — nl 2 


(2.89) 


According to p. 428 of Ref. [16], the angular velocity L u> B of B in L can be expressed as 

L u B = 0iS 2 + 0 3 )bi + (01 C 2 C 3 + e 2 S 3 ) b 2 + { 8 2 C 3 - 8 iC 2 S 3 )b 3 (2.90) 


Consequently, 

n uj b = ttibi + it 2 b 2 + u 3 b 3 

= JV w t + V (2.91) 

= ( 8^2 + e 3 - nS 2 ) bi + { 8 x 0203 + 8 2 S 3 - nC 2 C 3 ) b 2 + ( 8 2 C 3 - 8 iC 2 5 3 + nC 2 5 3 )b 3 

from which three kinematical equations can be obtained, 

a {C 3 U 2 ~ S 3 u 3 ) h a < r< a , S 2 {-C 3 u 2 + S 3 u 3 ) 

8 1 — h n, 8 2 — S 3 u 2 + C 3 u 3 , 8 3 — u\ H — (2.92) 

0 2 u 2 

If b r are taken to be parallel to central principal axes of inertia of S, as they are in the latter 


part of Sec. 2.2, then Eqs. (2.87) give rise to three dynamical equations, 

hui = (h - 7 2 )(3n 2 r 2 r 3 - u 2 u 3 ) - t x - f\ + w x (2.93) 

I 2 U 2 = (I\ — h){3n 2 rir 3 — u\u 3 ) — r 2 — r 2 + w 2 (2.94) 

hu 3 = {h - /i)(3n 2 rir 2 - u\u 2 ) — t 3 - t 3 + w 3 (2.95) 

where n = t • bj, n = r • bj {i = 1, 2, 3), and in view of Eqs. (2.88), 

ri = SiC 2 , r 2 = -(5i5 2 C 3 + 5 3 Ci), r 3 = 5i5 2 5 3 - C3C1 (2.96) 

The purpose of developing a control law will be to determine values of r* and fj (i = 1,2,3) that 
best control the orientation of B in L, and minimize the magnitudes of h and H. In addition 
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to Eqs. (2.92) and (2.93)-(2.95), six more equations will be used, three of which are simply a 
rearrangement of Eqs. (2.84)-(2.86), 


Hi = u 3 H 2 - U2H3 + fi ( 2 . 97 ) 

H2 = u\ H3 — U3H1 + T2 ( 2 . 98 ) 

H3 = U2H1 — U1H2 + T3 ( 2 . 99 ) 

and three of which are analagous and obtained from Eq. (2.75), 

hi = U3I12 - u 2 h3 + ti (2.100) 

h2 = Uih3 — U3hl + T2 (2.101) 

7-3 = U2hi — uih2 + T3 (2.102) 


2.5 Linearized Equations for Control Law Design 

A law for controlling the attitude of B with flywheels and CMGs can be developed from linearized 
forms of Eqs. (2.92)-(2.95) and (2.97)-(2.102). In the absence of w, t, and r (meaning, for the 
last two vectors, that the resultant of momenta h and H vanish, rotor spin speeds do not change, 
and CMG gimbals do not move) it is possible for B to remain fixed in L with l r having the 
same directions as b r , parallel to central principal axes of inertia of S'; one may linearize “about” 
this condition by letting the orientation angles 9 r (r = 1,2,3) have nominal values of zero, and 
introducing perturbations 9 1, 92 , and #3 such that 

0i = O + 6h, 0 2 = O + 0 2 , 03 = 0 + 03 (2.103) 

In addition, ui, U 2 , and U 3 are introduced as perturbations of the angular speeds of B in N, 

ui = 0 + ui, U 2 = -n + U 2 , U 3 = 0 + U 3 (2.104) 

Similarly, we linearize “about” H = h = 0 by introducing the following perturbations to H r and 
h r (r = 1,2,3), 

Hi = 0 + i?i, H 2 = 0 + H 2 , 773 = 0 + 773 (2.105) 

and 

7-1 = 0 + hi, 7-2 = 0 + 7-2, 73 = 0 + 7-3 (2.106) 

Differentiating Eqs. (2.103)-(2.106) with respect to time gives 

9 r = 9 r , u r = u r , H r = H r , h r = h r , (r = 1,2,3) (2.107) 

Kinematical differential Eqs. (2.92) are linearized by substituting from Eqs. (2.103) and (2.104), 
and discarding terms of second or higher degree in the perturbations. When the trigonometric 
functions of the attitude angles are expanded in power series, the linearization means 

sin 9 r = sin 9 r ~ 9 r , cos 9 r = cos 9 r « 1, (r = 1, 2, 3) (2.108) 
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Therefore, Eqs. (2.92) become 


h = [lfe - n) - 4%] + n 


« h 2 (2.109) 

#2 = # 3(^2 - n) + 1 % 

K,u 3 -n6 3 (2.110) 

~ _ 6 * 2 [-l(* 2 -n) + 0 3 u 3 ] 

C73 = H J 

«fti+n 0 2 ( 2 . 111 ) 

Similarly, linearization of Eqs. (2.93)-(2.95) performed with the aid of Eqs. (2.108), (2.104), 
and (2.96) yields 

hui = (h ~ h){3n 2 [6id 2 (l) + 0 3 ( 1 )][( 1 )( 1 ) - 9i0 2 9 3 ] - (u 2 - ra)(u 3 )} - n - n + uq 

~ ( h ~ h)(3n 2 9 3 + nu 3 ) - T\ - f\ + w\ (2.112) 

hu 2 = (I\ ~ T 3 ) {3tt, 2 (9i ( 1) [6*i 6* 2 (9 3 - (1)(1)] - uiu 3 } - r 2 - t 2 + w 2 

~ (h - /i)(3n 2 0i) - r 2 - f 2 + w 2 (2.113) 

hu:i = {I 2 - Ii){-3n 2 6'i(l)[0i6' 2 (l) + 0 3 (1)] - ui(u 2 - n)} -t 3 -t 3 + w 3 

~ {h ~ h)nui - t 3 - f 3 + w 3 (2.114) 

Upon appealing to Eqs. (2.104)-(2.107), the linearized forms of Eqs. (2.97)-(2.102) are 


H 1 = u 3 H -2 — («2 — n)ih 3 + n 
~ ni7 3 + n 

(2.115) 

H 2 = u\H 3 - u 3 H\ + r 2 
~ t 2 

(2.116) 

H 3 = (u 2 — n)Hi — u 1 H 2 + r 3 
~ —nH\ + r 3 

(2.117) 

hi = tt 3 h 2 - (u 2 - n)h 3 + T\ 
~ nh 3 + ri 

(2.118) 

h 2 = u\h 3 - u 3 h\ + r 2 
~ r 2 

(2.119) 

h 3 = ( u 2 - n)hi - uih 2 + r 3 
~ -nh\ + r 3 

(2.120) 


2.6 Nondimensional Equations for Control Law Design 


The linearized relationships (2.109)-(2.120) in Sec. 2.5 become better suited for control law design 
after they have been made nondimensional by use of the following definitions, 


0 * = 9i, 


* A Ui 

Ui = — , 
n 


, * A 

h* = 


I i n ’ 


H* = ^- (i = 1,2, 3) 
I 1 n 


(2.121) 
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* A Tj 
Ta — 


Ln 2 ' 


_* A Ti 
Ta — 


L n 2 ' 


* ^ ( • -| Q Q\ 

W i ~ t 2 (* — 1)2,3) 

In It 


.* A 

t — nt 


from which one obtains 


• a dd* ddi dt 9i 

= —*- = —±— = — (z = 1,2,3) 
dt* dt dt* n 

■ A du* d ( Ui\ dt m 

v* = - r ± = —[ — ]- = -£■ (z = 1,2,3) 

dt\n J dt* n 2 


dt* 


A 

dh* 

d( 

hi \ 

dt 

hi 


dt* 

dt y 

I,n) 

1 

1 * 
-+o 

1^3 

Un 2 

A 

dH* 

d( 


\ dt 

Hi 


dt* 

dt 1 


) dt* 

Un- 


(* = 1,2,3) 

(* = 1,2,3) 

Substitution from Eqs. (2.121), (2.122), and (2.124)-(2.127) into (2.109)-(2.120) yields 


and 


where ki in Eqs. (2.129)— (2.131) are defined as 


fci = 


A h — h 


h 


ko = 


a h — h 


h 


. a h — h 


(2.122) 

(2.123) 

(2.124) 

(2.125) 

(2.126) 
(2.127) 


= U * 2 , e* 2 = ut-et, o* = ui + o* 2 

(2.128) 

u\ = 3 fci 03 + k\U 3 — Ti — f{ + wl 

(2.129) 

u 2 = 3fc 2 0* — t 2 — f 2 + w 2 

(2.130) 

u* 3 = k 3 u{ - t 3 - f 3 + w 3 

(2.131) 

= K + Tt, h* 2 = t 2 , h* = -hl + Ti 

(2.132) 

+ ft, H* = f 2 , H* = -HI + f* 

(2.133) 


(2.134) 


Additional differential equations governing f hi dt and f Hi dt, such as those introduced in 
Ref. [6], are sometimes needed to eliminate biases in the time histories of hi and Hi (i = 1,2,3). If 
one makes the following definitions 


x i+ i 2 = J h* dt * , x i+15 = J H* dt* (i = 1,2, 3) 

then the associated differential equations are given by 


(2.135) 


x i+12 = dx i+12 /dt* = h*, Xi +15 = dx i+15 /dt* = H* (i = 1,2,3) 


(2.136) 
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If the control law is to maintain a specified orientation of B in L, feedback of integrals of attitude 
errors is required, therefore we define the quantities 

Xi+is = f 6*df (z = 1,2,3) (2.137) 

which are governed by the differential equations 

x i+18 = dxi+w/dt* = 6* (i = 1,2, 3) (2.138) 

The foregoing nondimensional, linear differential equations have the form required for design of 
linear control laws 

{x} = [A]{x} + [B]{r} + {W} (2.139) 

where {x} is a 21 x 1 column matrix whose first 12 elements are defined as 

Xi = e*, X i+3 = u*, x i+6 = h*, x i+9 = H*, (* = 1,2,3) (2.140) 

and whose final 9 elements are defined in Eqs. (2.135) and (2.137), respectively. The elements of 
the 21 x 1 column matrix {£} are defined as 

±i = dxi/dt* (i= 1, . . . , 21) (2.141) 

and the column matrix {r} is dimensioned 6x1 with the elements 

M = [ t* r 2 * t 3 * t* t| r 3 * ] T (2.142) 


where the superscript T indicates the transpose of a matrix. The 21 x 1 column matrix {IT} is 
defined as 


{IT} 

To 

0 

0 

wf 

w* 2 

w l 

0 

0 

0 

0 0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

] T 

(2.143) 

The 21 x 21 square matrix [.A] 

is 

















' 0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 ' 



0 

0 

-1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 



0 

1 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 



0 

0 

3Aq 

0 

0 

h 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 



3 k 2 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 



0 

0 

0 

k'3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 



0 

0 

0 

0 

0 
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0 
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0 

0 

0 

0 

0 

0 
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0 

0 

0 

0 
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0 

0 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 



0 

0 

0 

0 

0 

0 

-1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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[A] = 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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and the matrix [B] is dimensioned 21 x 6, 


[B] 


A 


0 

0 

0 

-1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0-10 0 
-10 0-10 
0-10 0-1 
0 0 0 0 0 

1 0 0 0 0 

0 10 0 0 

0 0 10 0 

0 0 0 1 0 

0 0 0 0 1 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 


(2.145) 


2.7 Nonlinear Equations used in Simulations 

In the physical system, relative motion between a flywheel rotor and the spacecraft will undoubtedly 
meet with some form of resistance that will necessitate consumption of power to overcome and, if 
not countered, will cause the actual amount of stored energy to diverge from the required amount. 
Section 4.3 deals with a control law designed to compensate for any such error in flywheel rotational 
kinetic energy, and consequently we now introduce terms to represent damping in the equations 
used for simulations of nonlinear motion, the results of which are presented in Chapters 4 and 5. 

Up to this point the moment about R* of forces exerted by B on Ri has been represented 
by M. B / Ri ] henceforth we regard the dot product M. B ^ Rl • j3 { (i = 1 , . . . ,T) as the sum of two 
contributions. The first is from a motor-generator, and will be denoted by • /3j for conve- 

nience. The second is due to damping, related by a constant of proportionality C d to the angular 
speed of Ri relative to B. For example, in connection with rotor i? 4 , • bq is replaced by 

M B / Ri • bi — CdU 5 in Eq. (2.55). As a result Eqs. (2.87) give way, after decoupling, to 

(I\ - 2J)u\ = (I 3 - / 2 )(3n 2 r 2 r3 - u 2 u 3 ) - J[u 2 {u 8 + ug) - ^ 3(^6 + u 7 )] 

+ W\ — 7~l — ( M b / r ' + M s / fi4 ) • bi + C d {u 4 + u 5 ) 

{h ~ 2 J)u 2 = (h ~ / 3 )(3n 2 nr 3 - uiu 3 ) - J[u 3 (u 4 + u 5 ) - ui(u 8 + «g)] 

+ U> 2 - T 2 - ( M b / R2 + M b / r • b 2 + C d (u 6 + u 7 ) 

(I 3 - 2 J)ii 3 = (I 2 - Ji)(3n 2 rir 2 - mu 2 ) - J[ui(u 6 + u 7 ) - u 2 (u 4 + u 5 )} 

+ w 3 -t 3 - (m b / r * + M b / r ^ ■ b 3 + C d (u 8 + ug) 


(2.146) 

(2.147) 

(2.148) 
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and Eqs. (2.54)-(2.59) governing flywheel speeds are restated as 


J(U! + u 4 ) = M B/Rl • bi - C d u 4 (2.149) 

J{u\ + U 5 ) = • bi — C d u§ (2.150) 

J(ii 2 + ue) = M b ! R2 • b 2 - CdUQ (2.151) 

J{u 2 + * 7 ) = M B / R5 • b 2 - C d u 7 (2.152) 

J(h 3 + u 8 ) = M b/Rs • b 3 - C d u 8 (2.153) 

J(u 3 + it 9 ) = M B / R * -b 3 -C d u 9 (2.154) 


where the dot products M,; • bj and Mj +3 • bj (i = 1, 2, 3) are neglected for the reasons discussed at 
the conclusion of Sec. 2.2. The simulation results presented in Chapters 4 and 5 are obtained from 
numerical solutions of nonlinear differential equations (2.146)-(2.154), together with the differential 
equations (2.100)-(2.102) for CMG momentum, kinematical differential equations for quaternion 
elements [see Eqs. (5.7) and (5.8)] in place of Eqs. (2.92), and dimensionalized counterparts to 
Eqs. (2.136) and (2.138) for integrals of momentum and attitude. 
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Chapter 3 

Flywheel Steering Laws 


The great benefit of utilizing flywheels is that they can serve simultaneously as attitude control 
actuators and as energy storage devices; this dual role requires that t, obtained on the basis of 
attitude control considerations, be applied in a way that allows energy to be stored or discharged 
as needed. A flywheel rotor Ri is suspended in a vacuum housing in B with magnetic bearings, 
and relative motion between Rj and B is brought about by a motor-generator that enables B to 
exert on Ri a torque with magnitude • /3j, the purpose of which is to produce and change 

rotor momentum in order to furnish attitude control, and to alter the rotor’s rotational kinetic 
energy. An important measure of energy storage is power, or the rate at which rotational kinetic 
energy is changed. In this chapter relationships for M B ^ R ' • /3j as functions of power, t, and 
rotor speeds are developed; they are referred to collectively as a flywheel steering law because they 
are similar in nature to a CMG steering law that determines girnbal speeds (and thus, indirectly, 
gimbal motor torques) needed to produce r as requested by a control law. Two such steering laws 
are presented; the first is the result of simply prescribing the total power of the flywheel system, 
whereas specification of the power required for each of three flywheel pairs gives rise to the second 
law. Bearing friction and damping are neglected in the design of the steering laws. 

The flywheel rotors are arranged in counter-rotating pairs; each pair, denoted by E), consists of 
rotors Ri and Ri + 3 (i = 1,2,3). Referring to Eq. (5) of Ref. [22], the power B P Fi of T) in B can 
be expressed as 

B P Fi = (i R i/ R t . B a Ri ^j ■ B ua Ri + (l R ’+ 3 / R h3 . B a Ri + 3 ^j ■ B u Ri + 3 (i = 1,2,3) (3.1) 

or, in view of Eqs. (2.36) and (2.37), 

B P F " = J (u2i + 2U2i+2 + U2i+3ll2i+3 ) (* = 1 , 2 , 3) (3.2) 

and therefore the total power may be expressed as 

3 3 

B p t B P Fi = ( U 2 i+ 2 U 2 i +2 + U2i+3U2i+3) (3.3) 

2=1 2=1 

Now, if iii, ii 2 , and 113 are assumed to be small in comparison to 114 , . . . ,uq, and if M, • bj 
and M i+3 ■ b* (i = 1,2,3) vanish for the reasons put forth at the conclusion of Sec. 2.2, then 
Eqs. (2.54)-(2.59) may be approximated as 

Jii2i+2 ~ M B/Rl • bj, Jii2i+3 ~ M B/Ri + 3 • b* (z = 1 , 2 , 3) (3.4) 
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which, together with Eqs. (2.81) and (2.83), lead to 


(r - n ld b x H) • b i = Hi = J(u 2 i +2 + ^ 24 + 3 ) 


(M B/Ri + M B/Ri+3 ^j • bi (i= 1,2,3) 


or 


M B/R l+3 . bi » - N UJ B x H - M B / Ri J • b i (i = 1,2,3) 

Substitution from Eqs. (3.4) into Eq. (3.3) produces 


(3.5) 


(3.6) 


B ; 


~ E («2* 2 M*/* + U2i+3 M s /^+3) • b 


(3.7) 


i— 1 


Eqs. (3.6) and (3.7) constitute a system of four equations, linear in the six unknowns M B/ R ' • b, 
and ^sA B / Ri + 3 • b^, (i = 1,2,3). This underdetermined system can be written in matrix form as 


where 


[M { y } = (4 



1 1 0 0 0 0 
0 0 110 0 
0 0 0 0 1 1 




144 U-5 

uq 

u 7 

14 S U9 

m b / Ri 

• bi 

' 



M-B/i? 4 

• bi 




M B / R2 

• b 2 





• b 2 


> 


M B / Rs 

• b 3 




M B /«6 

• b 3 




(r- N 

u; B 

X 

H) 

• bi 

e f- N 

us B 

X 

H) 

• b 2 

(r-- 

i o B 

X 

H) 

- 

• b 3 

B P 

F 



(3.8) 


(3.9) 


(3.10) 


(3.11) 


One may solve Eqs. (3.8) by forming a matrix pseudo-inverse such as the one presented in 
Ref. [23] and developed by Moore and Penrose for underdetermined systems, 

[^] + = [^] T ([^][^] T ) _1 (3.12) 

which yields the solution 

{#} = [4J + {4 (3.13) 

that minimizes the sum of the squares of the unknowns, {y} T {y}. (There exist several other 
performance measures that could be considered in solving an underdetermined system of equations.) 
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Use of a pseudo- inverse is in essence the suggestion made by Hall in Sec. 4 of Ref. [9], yielding a 
steering law wherein the power and attitude control requirements are met simultaneously, and a 
function of the instantaneous motor torques, the sum • bj) 2 + (M iJ// ~ Ri+3 • b i) 2 ], is 

minimized. 

Upon making the definitions 


A 

= U5- U 4 , 

, A 

d 2 = U7 - UG, 

00 

llt> 

A 

= U 5 + U 4 , 

A 

s 2 = U 7 + U 6 , 

A 
S3 = 


(3.14) 

(3.15) 


the pseudo-inverse can be written explicitly as 


[A a ] + = (3.16) 



2rt5<ii + c?2 2 + d 3 2 

s 2 di 

s 3 di 

— 2di ' 


— 2u4<ii + c?2 2 4" d 3 2 

-s 2 di 

-s 3 di 

2di 

1 

sid 2 

2u 7 d2 + di 2 + d 3 2 

s 3 d 2 

—2 d 2 

2 (d! 2 + d 2 2 + ds 2 ) 

~s\d 2 

—2ugd2 + di 2 + d3 2 

-s 3 d 2 

2 d 2 


sid 3 

S2d 3 

2ugd 3 + di 2 + d 2 2 

-2 d 3 


-sid 3 

~S2d 3 

—2usd3 + d 3 2 + d 2 2 

2 d 3 J 


It is worth noting that if the rotor speed differences d{. o ?2 , and d 3 vanish for all three flywheel 
pairs, the pseudo- inverse becomes infinite and the steering law does not furnish a result. Since the 
rotor speeds of the flywheels in each pair will normally have opposite signs, this condition should 
be unlikely. In addition, examination of Eqs. (3.10), (3.11), (3.13), and (3.16) reveals that when 
the flywheels are not required to provide attitude control (r = 0), and the condition of counter- 
rotation is present ( u 2 j +3 = — U 21+2 , thus H = 0), this steering law dictates that Wl B / Ri + 3 • bj = 
— M B / Ri • bj, and counter-rotation is preserved. 

As an alternative to dealing with the four Eqs. (3.6) and (3.7), one could replace the single 
Eq. (3.7) with three others by choosing to divide the power requirement evenly among the three 
flywheel pairs, B P 1 ' 1 = B P F ' 2 = B p }, 3 — iBpF> yi e pp n g s j x equations in six unknowns. Substi- 
tution from Eqs. (3.4) into Eqs. (3.2) gives the three new equations 


B P Fi « (u 2 i+- 2 M B,Ri + U 2 i+m B,Rt+3 ) ■ bj (i = 1,2, 3) (3.17) 


Now, substitution from Eqs. (3.6) gives 


B P Fi « U2i+2 ( M B / Ri • bj) + U 2 i+3 (t — N (jD B X H - M b / r <) • b t 

= ( U2i+2 ~ U2i + z) (m b/r ‘ • bj) + U2i+3 [j ~ N ^J B X H) • bj (i = 1, 2, 3) (3.18) 

which can be rearranged to yield 


M B/Ri ■ bj = 


B P Fi - U 2 i+3 (t - N, ' B 


X 


H) - bj 


U2i+2 ~ U2i+3 

and, when one substitutes from this expression into Eq. (3.6), the result is 

H) • bj - B P Fi 


(* = 1,2,3) 


M B/R i+3 . ^ . = 


U2i+2 (T - N U B X 


U2i+2 ~ U2i+3 


(* = 1,2,3) 


(3.19) 


(3.20) 
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Eqs. ( 3 . 19 ) and ( 3 . 20 ) constitute another flywheel steering law, indicating the moment that must 
be applied by a motor-generator to each of two rotors belonging to a counter-rotating pair in order 
to apply t as called for by an attitude control law, and at the same time satisfy power requirements 
specified by B P Fi . Each flywheel pair is expected to operate with the sign of the rotor speed U2i+3 
opposite the sign of U21+2 hence, the denominators U2i+2 ~ U2i+3 should remain well away from 
zero. In the event that the flywheels are not required to participate in attitude control, and the 
condition of counter-rotation is present, the steering law yields ^A B / Ri + 3 • bj = — ^J[ B / Ri . and 
thus maintains the condition of counter-rotation. 
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Chapter 4 

Linear Controllers 


The linearized, nondimensional equations presented in Sec. 2.6 are used together with the Linear 
Quadratic Regulator (LQR) technique to obtain linear laws for controlling motion associated with 
two important orientations of an Earth-pointing spacecraft: torque equilibrium attitude (TEA), 
and attitude that is fixed (or held) with respect to a local-vertical-local-horizontal reference frame. 
Sections 4.1 and 4.2 contain the numerical values for LQR weighting and other simulation parame- 
ters, respectively. The design of a method for feeding back error in the rotational kinetic energy of 
the flywheel rotors in order to eliminate problems caused by rotor damping is discussed in Sec. 4.3. 
The comprehensive control system formed by the combination of a LQR for attitude control and 
momentum management, a flywheel steering law, and feedback of kinetic energy error is depicted 
in Sec. 4.4 with a block diagram. Section 4.5 begins with a brief discussion of the design of a control 
law that seeks TEA with CMGs and flywheels, and continues with a presentation of simulation 
results that illustrate the performance and features of the control and steering laws. First, TEA 
seeking behavior in the absence of flywheel rotor damping is examined, then the consequences of 
failing to counteract damping, and the advantages of kinetic energy error feedback are demon- 
strated in turn, and finally the performances of the pseudo-inverse and divided power steering laws 
are compared in some detail. The subject of Sec. 4.6 is the design and application of several control 
laws for holding attitude: without management of angular momentum of flywheels or CMGs, with 
inclusion of flywheel momentum management, and by means of thrusters. In conclusion, a com- 
parison is undertaken of the amount of propellant required if each of these operational techniques 
is used in connection with the rendezvous and docking of ISS with another spacecraft. 

4.1 Weighting 

Laws for controlling the orientation of B in L, and the momentum of flywheels and CMGs, have 
been designed using the infinite-horizon Linear Quadratic Regulator technique, in which a scalar 
quadratic performance index given by 

V = jf” ({x} T [Q] {*} + {t} t [B] {t}) dt (4.1) 

is minimized subject to the linear Eqs. (2.139) with {W} = {0}. The technique yields a state 
feedback gain matrix [K] which in turn is used to obtain {r}, 

M = - [K] {x} (4.2) 
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Table 4.1: State Weightings 


Parameter 

Max Value 

Nondimensional Value 

i = 

1,2,3 

1 

2 

3 

0i 

1 deg 

3282.8 

3282.8 

3282.8 

Ui 

0.2 deg/s 

0.1050 

0.1050 

0.1050 

hi 

6,779 N-m-s 

237.9 

11.0 

322.9 

Hi 

6,779 N-m-s 

237.9 

11.0 

322.9 

f hidt 

2.7 x 10 6 N-m-s 2 

1161.7 

53.6 

1576.9 

f Hidt 

2.7 x 10 6 N-m-s 2 

1161.7 

53.6 

1576.9 

J 6idt 

1000 deg-s 

2564.7 

2564.7 

2564.7 


As suggested by Bryson and Ho in Ref. [24] , the weighting matrix \Q\ can be chosen as diagonal, 
and unity should be approximately equal to the product of Qjj and the square of the maximum 
acceptable value of the associated element Xj of {x}, (j = 1, . . . ,21). Likewise, a diagonal form 
of [R] is convenient, with unity approximately equal to the product of Rkk and the square of the 
maximum acceptable value of the associated element of {r} (k = 1 , . . . ,6 ). 

Maximum acceptable values of several parameters used in constructing \Q\ are listed in Table 
4.1. A maximum maneuver rate of 0.1 deg/s specified on p. 7-6 of Ref. [25] provides a limit on the 
magnitude of L u) B \ since the orbital angular speed n will not be greater than about 0.065 deg/s, 
a limit of 0.2 deg/s imposed on each of Hi appears reasonable. Each CMG possesses 4,745 N-rn-s 
(3,500 ft-lbf-s) of momentum, therefore the magnitude of the resultant momentum of 4 CMGs is 
limited to 18,981 N-rn-s (14,000 ft-lbf-s). Even if each of three components of momentum reach an 
absolute value of 6,779 N-rn-s, the magnitude is well within the limit. For the sake of simplicity, 
the same limits are applied to flywheel momentum; it should be kept in mind that the spin speeds 
of the International Space Station flywheels are limited in absolute value to 60,000 rpm, the axial 
central principal moment of inertia J is expected to be approximately 0.3010 kg-m 2 , and there 
will be no more than 96 flywheel rotors (arranged in 48 pairs). Limits of 2.7xl0 6 N-m-s 2 on the 
integrals of momentum correspond to keeping momentum in each axis below 544 N-m-s for 5,000 
s (almost one orbit). Similarly, the limit of 1000 deg-s on the integrals of attitude error arises 
from requiring the attitude error to be less than 0.2 deg for 5,000 s. Before constructing [Q], the 
values in the second column of Table 4.1 are made nondimensional to be in correspondence with 
the elements of {x}, as indicated in the last three columns. 

The maximum expected value of each element of {r} is taken to be 2.7 N-m; nondimensional- 
ization yields Rn = R 44 = 1903.2, R 22 = R 55 = 87.8, and R33 = Rqq = 2583.6. 

It is important to point out that the preliminary control law designs presented herein do not 
address performance and stability robustness in the presence of model uncertainty. Such detailed 
designs are beyond the scope of the present work, although they would be required for an actual 
spacecraft. 
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4.2 Simulation Parameters 


Other values required for a numerical simulation are as follows. 

Moments and products of inertia of S with respect to S* are taken to be 




50.28 

-0.39 

-0.24 ' 

'jS/S *' 

= 1.36 x 

-0.39 

10.80 

0.16 



-0.24 

0.16 

58.57 


x 10 6 kg-rn 2 


(4.3) 


which are the same values (in metric units) as those associated with Phase 1 in Table 1 of Ref. [6] , 
with the exception of /13 which is suspected to be a typographical error, the correct value being 
—0.24 x 10 6 slug-ft 2 . 

The value of w used in Ref. [6] represents the moment about S* of aerodynamic forces exerted 
on the Phase 1 configuration of S when the attitude is near TEA; it is given by 

w = 1.36[(1 -f- sin nt H — sin2nf)bi + (4 + 2 sin nt + - sin2nf)b>2 + (1 + sin nt H — sin2nt)b3] N-m 

(4.4) 


where n, the magnitude of N u L , is taken to be 0.001131 rad/s. 

Each flywheel pair in the physical system is required to discharge 4,400 W of power during 
the portion of the orbit that lies within the Earth’s shadow, known as the period of eclipse. The 
remaining portion of the orbit, during which sunlight reaches the spacecraft, is taken to be twice 
as long as the eclipse. Therefore, for each orbit, the total power that must be supplied by the 48 
pairs of flywheels in the physical system is given by 


BpF = 


48 x 2, 200 W = 105.6 kW 

48 x -4, 400 W = -211.2 kW 


0<t<~ (charge) 

(discharge) 


(4.5) 


where the bar over P indicates a known function of t to be used in connection with the pseudo- 
inverse steering law developed in Ch. 3. The alternative steering law is referred to as a divided 

power law because the power requirement is divided into three equal parts, therefore B P hl = 
BpF 2 = BpF 3 = 1 BpF or 


35.2 kW 0 <t< (charge) 

— 70.4 kW < t < p- (discharge) 


(4.6) 


As mentioned earlier, the physical system is made up of 96 flywheels; since the present model 
involves only 6 rotors, a scaling factor of 96/6 = 16 is used; thus J = 16 x 0.3010 kg-rn 2 = 4.82 
kg-rn 2 . 


4.3 Energy Feedback 

The power possessed by the flywheels, given in Eq. (3.3), can differ from the value required by 
Eq. (4.5) if, for example, a rotor Ri and B exert on each other a couple whose torque is proportional 
to the rotor speed. If left unchecked, such unwanted resistance or damping will lead to a difference 
between B p 1, and B P F that increases with time, to rotor speeds that exceed their maximum and 
minimum limits, and to singularities in the steering laws. These deleterious effects of damping can 
be eliminated by the control system through feedback of rotational kinetic energy error. 
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The total rotational kinetic energy B K F of the flywheel rotors relative to B can be expressed 
as 

b k f = ^Y, u ? ( 4 - 7 ) 

i — 4 

and the power of F in B is given by the derivative of B I\ F with respect to t, 

BpF = d BR F 

dt v ’ 

One can regard the required power B P F as the time derivative of a required kinetic energy of F 

in B 

B p F = — B K F (4.9) 

dt K ’ 

and define a kinetic energy error e k as the quantity 

e k = B K F - B K F (4.10) 

that is governed by the differential equation 

e k = — ( b K f - b K f ) = B P F - B P F = e„ (4.11) 

dt x 

where e p is defined to be the error in power, or the difference between the actual and required 
values. The LQR technique can be used to control the kinetic energy error by minimizing the cost 
function 

r oo 

Ve = / (A e k ~ + e p 2 ) dt (4-12) 

Jo 

where A is a weighting parameter on the kinetic energy error of the flywheel system. This leads to 
a feedback controller in which the required power is adjusted by the amount —y/\e k . Thus, the 
commanded power is defined to be 


B P c F = B P F -V Xe k (4.13) 

and is used in place of B P F in Eq. (3.11) for the pseudo-inverse steering law. Similarly, a com- 
manded power B P Fi (i = 1, 2, 3) is obtained for each of the three flywheel pairs and used together 
with Eqs. (3.19) and (3.20) for the divided power steering law. The merits of this kinetic energy 
error feedback are illustrated presently. 

4.4 Control System Block Diagram 

The block diagram in Fig. 4.1 shows the LQR (controller), flywheel steering law, and kinetic 
energy error feedback arranged to form a complete control system in MATLAB/Simulink®. The 
differential equations governing the behavior of the plant model (see Sec. 2.7) are placed in a 
Simulink S-Function to be numerically integrated and provide the value of the state {x} at the 
current simulation time. The current value of the state is fed to the controller S-Function, where 
it is made nondimensional and then multiplied by the LQR-derived state feedback gain matrix to 
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Figure 4.1: Control System Block Diagram 


generate the control torques r and t. The controller also determines the kinetic energy error of 
the flywheel system, e*,, needed to calculate the commanded power, B P C F , that is used together 
with the flywheel control torque f in the flywheel steering law to obtain the flywheel motor torque 
measure numbers in {y} [see Eq. (3.10)]. These are supplied to the plant model, and the control 
loop cycle is repeated. (The subject of CMG steering is well understood; therefore, our simulations 
do not include a model of individual CMG gimbal motions.) 

4.5 TEA Seeking 

The Torque Equilibrium Attitude (TEA) of a spacecraft is defined as the orientation for which 
the angular acceleration of B in N vanishes. The present ISS attitude and CMG momentum 
control algorithm can keep the orientation in the neighborhood of a time-varying TEA almost 
indefinitely, without requiring any expenditure of propellant from the Reaction Control System. In 
what follows, the design of a control law that seeks TEA with CMGs and flywheels is described 
briefly, and then simulation results are presented that illustrate the performance and features of 
our control and steering laws. We examine TEA seeking behavior in the absence of flywheel rotor 
damping in Sec. 4.5.1, illustrate the consequences of failing to counteract damping in Sec. 4.5.2, and 
then demonstrate the advantages of kinetic energy error feedback in Sec. 4.5.3. The performance 
of the pseudo-inverse steering law is compared to that of the divided power scheme in Sec. 4.5.4. 

Seeking a TEA requires regulating the states contained in the 18 x 1 column matrix 

M = 0* 2 , 6*, hi,..., HI,..., Jh*dt*,..., fH* 1 dt*,...] T (4.14) 

Average TEA angles are in general nonzero; therefore, the state does not include f Qidt as elements. 
Weighting matrices [Q] and [R] are constructed with the associated values presented in Sec. 4.1. 
The matrices [A], \Q\, and [B] are dimensioned 18 x 18, 18 x 18, and 18 x 6, respectively, in 
accordance with the number of regulated state variables for this control scheme. 

Figure 4.2 shows the eigenvalues of the matrix [A] — [ B ] [K], or the closed-loop poles of the 
system. The diagonal lines mark the points for which the magnitude of the eigenvalue’s real 
number is equal to the magnitude of the eigenvalue’s imaginary number; in other words, the points 
at which a damping ratio of l/\/2 is obtained. Regional pole-placement methods presented in 


37 


TEA-Seeking Closed-Loop Poles 



Figure 4.2: TEA-Seeking Closed-Loop Poles 


Ref. [26] may also be used to generate pole locations resulting in damping ratios greater than or 
equal to l/\/2; however, such methods are beyond the scope of the current research. The transpose 
of the 6 x 18 gain matrix [K] associated with these poles is given in Eq. (4.15). 


- 0.0000 

-4.8044 

- 0.0000 

- 0.0000 

-4.8044 

- 0.0000 

-1.8758 

- 0.0000 

-0.9269 

-1.8758 

- 0.0000 

-0.9269 
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(4.15) 


Simulation results for this control law, discussed in the remainder of this section, are obtained 
with the following initial values of the state variables. Angles describing the orientation of B in L at 
t = to are 6*i (to) = 5° (pitch), # 2 (^ 0 ) = 5° (yaw), and # 3 (to) = 5° (roll). Angular speeds associated 
with N u> B (with L u B = 0 ) are ■ui(fo) = —9.86 x 10~ 5 rad/s, ^ 2 (^ 0 ) = —1-12 x 10~ 3 rad/s, and 
us(to) = 9.82 x 10— 5 rad/s. Rotor spin speeds are « 4 (to) = U(i{to) = us(to) = —20,000 rpm, and 
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us(to) = «7(io) = ug(to) = 20,000 rpm. Initial values of CMG momentum measure numbers are 

hi(to) = h 2 {t 0 ) = hs(to) = 0. 

4.5.1 Undamped Flywheel Rotors 

Figures 4. 3-4. 8 show performance of the TEA-seeking control law with undamped flywheel rotors, 
obtained by setting Cd = 0 in Eqs. (2.146)-(2.154). 

A time history of the orientation of B in L, described with a body-three, 2-3-1 sequence, is 
shown in the upper plot of Fig. 4.3; the solid curve is used for 6 \ (pitch), the dashed curve for d 2 
(yaw), and the dash-dot curve for 63 (roll). The average values of these orientation angles in the 
steady state are referred to as average torque equilibrium attitude angles, and are approximately the 
same as those shown in Ref. [6], —7.5° in pitch, —1.2° in yaw, and —0.2° in roll. The amplitudes of 
the steady state oscillations can be reduced significantly with cyclic disturbance rejection filters, as 
shown in Ref. [6]. The initial values of the angles, being rather different from the torque equilibrium 
attitude, lead to the transient behavior evident in the first few orbits in Figs. 4. 3-4. 8. 

The lower plot of Fig. 4.3 illustrates the behavior of u\, u 2 , and 113 , shown with solid, dashed, 
and dash-dot curves respectively; their respective average values are 0, — n, and 0, the condition 
about which linearization was performed [see Eqs. (2.104)]. 

Figure 4.4 illustrates the progression of the magnitudes of CMG momentum h in the upper 
curve, and flywheel momentum H in the lower curve. The former is obtained from numerical 
integration of Eqs. (2.100)-(2.102), whereas the latter is calculated according to Eq. (2.39). In the 
latter portion of the simulation, CMG momentum magnitude remains well below 4,745 N-rn-s, the 
capacity of a single CMG. The similarity of the two curves in Fig. 4.4 is to be expected in view of 
the equal weights on their errors used to obtain the gains. 

Time histories of measure numbers of so-called CMG and flywheel control torques, t and r, are 
presented in Fig. 4.5; t\, t 2 , and 73 are shown with solid, dashed, and dash-dot curves respectively 
in the upper plot, whereas ?i, f 2 , and T3 appear in the lower plot. In the steady state, all of these 
quantities have reasonable amplitudes of 3 or 4 N-m; since the CMG and flywheel momenta are 
weighted equally, it is not surprising that the curves in the upper plot appear to be the same as 
the corresponding curves in the lower plot. 

Figure 4.6 displays the motor torques that must be applied to the flywheel rotors in order to 
meet the attitude control and power management requirements. Motor torques of the pair whose 
spin axes are parallel to bi are contained in the upper plots, with fA B ^ Rl • bq represented by the 
solid curve, and M Sy/i?4 • bi indicated by the dashed curve. Motor torques of the pairs whose spin 
axes are parallel to b2 and b3 are shown in the middle and lower plots, respectively. The first 
column of Fig. 4.6 shows the flywheel motor torques as a function of time using the pseudo-inverse 
steering law, the first of two methods presented in Ch. 3, and the second column shows the response 
using the second method, referred to as the divided power steering law. The time histories resulting 
from the two steering laws are comparable, and both ultimately lead to the same control torque t 
and required power B because this is the objective of both steering laws. The discontinuities 
in both instances occur as a result of alternating between charging and discharging in the required 
power profiles given by Eqs. (4.5) and (4.6). 

Figure 4.7 displays the angular speeds of the 6 flywheel rotors relative to B. Speeds of the pair 
whose spin axes are parallel to bi are contained in the upper plots, with 14 represented by the 
solid curve, and 74 indicated by the dashed curve. Speeds of the pairs whose spin axes are parallel 
to b2 and b3 are shown in the middle and lower plots, respectively. The first column of Fig. 4.7 
shows the flywheel angular speeds obtained with the pseudo-inverse steering method, and second 
column of Fig. 4.7 shows the response using the divided power steering law. The angular speeds 
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Figure 4.3: Attitude, Angular Velocity 
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Figure 4.4: Resultant Momentum Magnitude 
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Figure 4.5: Control Torques 
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Figure 4.6: Flywheel Motor Torques (Pseudo-Inverse left, Divided Power right) 
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.7: Flywheel Angular Speeds (Pseudo-Inverse left, Divided Power right) 
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Figure 4.8: Power, Power Error, and Kinetic Energy Error 
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are similar, but not identical, because the flywheel motor torques are not unique, even though both 
steering laws produce flywheel motor torques that will meet the attitude control torque and power 
management requirements. With the exception of us during a brief period in the first orbit, the 
absolute values of the rotor speeds remain less than 60,000 rpm, as they must. In keeping with 
typical flywheel designs, all of the rotor speeds have a ratio of maximum to minimum absolute 
value of approximately 3 to 1 as they charge and discharge. 

The total power of the flywheels B P F is calculated from the motor-generator torques and wheel 
speeds by means of Eqs. (3.7), and is seen in the upper plot of Fig. 4.8 to have very nearly the profile 
specified in Eqs. (4.5), 105.6 kW when the flywheels charge, and -211.2 kW when they discharge. 
In fact, the second and third plots of e p and e *. [see Eqs. (4.10) and (4.11)] reveal no errors in power 
or kinetic energy; this is to be expected since Cd = 0. It is clear that the power tracking objectives 
are achieved at the same time that the flywheels are used to assist the CMGs in attitude control. 

4.5.2 Damped Flywheel Rotors 

The detrimental effects of damping are brought to light by repeating the simulation described in 
Sec. 4.5.1, with an illustrative value of 10 “ 5 N-rn-s for the damping coefficient Cd- One might be 
tempted to neglect such a seemingly small effect, especially over the short term, but it is shown 
here to be troublesome if not dealt with over long periods. The performance of the TEA-seeking 
control law without feedback of kinetic energy error is recorded in Figs. 4.9-4.12. 

Figure 4.9 shows the time history of the attitude angles and the inertial angular velocity. As in 
the previous example, the solid curve is used for 6 \ (pitch), the dashed curve for 62 (yaw), and the 
dash-dot curve is used for 63 (roll). The attitude response is virtually the same, and the average 
TEA found by the control system is the same as in the case with undamped rotors. The lower 
plot of Fig. 4.9 shows the inertial angular velocity response, with u\, 112 , and 113 , shown with solid, 
dashed, and dash-dot curves respectively. 

The upper plot of Fig. 4.10 contains the actual power delivered to the spacecraft by the flywheel 
system as a function of time. The middle plot displays the error between the actual power and the 
required power due to the damping in the flywheel system, which leads to the large secular kinetic 
energy loss of more than 50,000 kJ after 10 orbits, shown in the lower plot. 

The kinetic energy error of the flywheel system is reflected in the angular speeds of the flywheels 
relative to B shown in Fig. 4.11. As in Fig. 4.7, speeds of the flywheel pair whose spin axes are 
parallel to tq are contained in the upper plots, with U 4 represented by the solid curve, and u§ 
indicated by the dashed; speeds of the pairs whose spin axes are parallel to t >2 and t >3 are depicted 
in the middle and lower plots, respectively. The flywheel rotor speeds exhibit a secular decay which 
is not present in the undamped flywheel speeds seen in Fig. 4.7. No significant difference appears 
between the behavior resulting from the pseudo-inverse and divided power steering laws, presented 
in the left and right columns, respectively. 

Time histories of flywheel motor torques, in the presence of damping, are displayed in Fig. 4.12 
where the results in the left and right columns are associated with the pseudo-inverse and divided 
power steering laws, respectively. A comparison with Fig. 4.6 reveals the effects of damping to be 
negligible at first; however, after some time it is apparent that damping causes the motor torque 
magnitudes to increase, with either steering law. Inspection of Eq. (3.7) indicates that the secular 
decay in rotor speeds requires an increase in motor torques in order to produce the required power. 
The increase in motor torque magnitudes implies further power losses as more electrical power must 
be diverted to the motor-generators in order to meet the attitude control and power management 
requirements simultaneously. 
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Figure 4.9: Attitude, Angular Velocity with Damped Flywheel Rotors 
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Figure 4.10: Power, Power Error, Kinetic Energy Error with Damped Flywheel Rotors 
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Figure 4.11: Flywheel Angular Speeds with Damped Flywheel Rotors (Pseudo-Inverse left, Divided 
Power right) 
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Figure 4.12: Flywheel Motor Torque with Damped Flywheel Rotors (Pseudo-Inverse left, Divided 
Power right) 
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4.5.3 Counteracting Damping with Kinetic Energy Error Feedback 

To compensate for damping, the kinetic energy error feedback designed in Sec. 4.3 is employed with 
a weighting parameter A of 1 s“ 2 . The parameters of the simulation discussed in Sec. 4.5.2 are used 
again, leading to the results reported in Figs. 4.13-4.15. 

Figure 4.13 shows that kinetic energy feedback eliminates the secular decay of rotor speeds seen 
to result from flywheel rotor damping in Fig. 4.11. 

The flywheel motor torques are shown in Fig. 4.14. It is immediately clear that the magnitudes 
do not increase with time as they do in Fig. 4.12. 

The power error e p shown in the middle plot of Fig. 4.15 is quite small and leads to the kinetic 
energy error displayed in the lower plot, which is small and periodic, in contrast to the secular 
decay obtained without energy feedback (A = 0 s~ 2 ). 

As evidenced by the results shown in this section, the kinetic energy error feedback method 
compensates for damping very effectively, and leads to power management results that are virtually 
the same as those for the undamped rotors. 


Flywheel Rotor Speeds 


Flywheel Rotor Speeds 







0 2 4 6 8 10 

Time [orbits] 


4 6 8 10 

Time [orbits] 


Figure 4.13: Flywheel Angular Speeds with Damped Flywheel Rotors and Kinetic Energy Error 
Feedback (Pseudo-Inverse left, Divided Power right) 


46 




Time [orbits] 


20 
10 
o 

-10 
-20 

0 2 4 6 8 10 

Time [orbits] 



Figure 4.14: Flywheel Motor Torques with Damped Flywheel Rotors and Kinetic Energy Error 
Feedback (Pseudo-Inverse left, Divided Power right) 
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Figure 4.15: Power, Power Error, and Kinetic Energy Error with Damped Flywheel Rotors and 
Kinetic Energy Error Feedback 
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4.5.4 Comparison of Steering Laws 

In the results presented thus far, two steering laws have been used to determine the flywheel motor 
torques necessary to meet the attitude control and power management requirements simultaneously. 
The first of these involves formation of a pseudo-inverse to solve an underdetermined system of 
equations for flywheel motor torques, whereas the second divides the power requirements evenly 
among the three pairs of rotors, resulting in a uniquely determined solution. The objective of both 
approaches is to produce the same control torque r and required power B P F ; therefore, it is not 
surprising that large differences are not discernible in the left and right columns of the figures 
presented in Secs. 4. 5. 1-4. 5. 3. A more detailed comparison of the differences is now undertaken 
with the aid of Figs. 4.16-4.18, where results obtained with the pseudo-inverse and divided power 
methods are displayed with solid and dashed lines respectively. 

Figure 4.16 shows time histories of the square root of the sum of the squares of the motor 
torques, {X]i=i • bj ) 2 + (M . B ^ Ri + 3 • bj) 2 ]}2, a quantity that is minimized by the pseudo- 

inverse approach. Curves associated with undamped flywheel rotors are contained in the first row, 
damped rotors in the second row, and damped rotors with kinetic energy error feedback in the third 
row. The solutions of the pseudo-inverse and divided power methods are indistinguishable from one 
another when viewed over an interval of 10 orbits in the left column; consequently, detailed plots 
are shown in the right column for the first half of an orbit. As expected, it can be seen that the 
minimum-norm solution of the pseudo-inverse method is always less than or equal to the result given 
by the divided-power method. These plots also seem to indicate that the pseudo-inverse solution 
approaches the divided-power solution in the steady state. The benefit of the kinetic energy error 
feedback is apparent once more, by comparing the growth in the second row of the first column 
with the absence of growth in the third row. 


Magnitude of Motor Torque 



8 10 




Magnitude of Motor Torque [Detail] 





Figure 4.16: Magnitude of Flywheel Motor Torque 
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Figure 4.17: Power Distribution 
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Figure 4.18: Energy Distribution 
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Figure 4.17 shows the distribution of power among the flywheel pairs, in the absence of damping. 
The first row contains the power B P Fl [see Eqs. (3.17)] of the rotor pair whose spin axes are aligned 
with bi, and the second and third rows show B p t2 and ^P^ 3 . A detail of the plots for the first 
2 orbits is shown in the right column of Fig. 4.17. 

Similarly, Fig. 4.18 displays the rotational kinetic energy possessed by each flywheel pair, 

in the absence of damping. The first, second, and third rows contain, respectively, B K Fi = 
J\(v, 2 i+ 2 ) 2 + ('W2i+3) 2 ]/2, (i = 1,2,3). The kinetic energies resulting from the two steering laws are 
indistinguishable when viewed over 10 orbits in the left column; therefore, detailed plots are shown 
in the right column for a single orbit. 

Naturally, the divided-power steering law distributes the power and energy equally between the 
three flywheel pairs. The pseudo-inverse steering law apportions the power profile among the three 
pairs in such a way that the total power requirement is met, but the minimum- norm solution for 
motor torques yields an unequal distribution of power and energy. As noted in connection with 
Fig. 4.16, the power and energy distributions resulting from the pseudo-inverse steering law appear 
to approach those of the divided power law in the steady state. Similar comparisons involving 
damping with and without kinetic energy feedback support these conclusions. Failure to meet 
power requirements exactly is remedied by the feedback of kinetic energy error, as previously 
discussed. 

4.6 Attitude Hold 

The operation involving rendezvous and docking of ISS with another spacecraft requires the at- 
titude of the ISS to be held fixed in a local-vertical-local-horizontal reference frame, an objective 
that differs from maintenance of a continuously varying TEA, and management of angular mo- 
mentum. Rather than continuously varying the spacecraft attitude to manage angular momentum, 
the controllers can be configured to keep the spacecraft in a specified orientation. In this section, 
several control laws for holding attitude are designed, and applied in simulations. Attitude hold 
without management of angular momentum of flywheels or CMGs is the subject of Sec. 4.6.1, the 
addition of flywheel momentum management is taken up in Sec. 4.6.2, and attitude hold by means 
of thrusters is discussed in Sec. 4.6.3. Finally, a comparison of the amount of propellant required 
by various operational techniques is presented in Sec. 4.6.4. In the interest of simplicity, only the 
pseudo-inverse steering law is exercised in performing these simulations, and flywheel rotor damping 
is not present. 

4.6.1 Attitude Hold without Momentum Management 

The LQR theory is applied here to the problem of holding a specified attitude without regard for 
the angular momentum possessed by flywheels or CMGs. The regulated states are 

{x} = [et, e* 2 , e* 3 , ut, u* 2 , u* 3 , feidt*, Jd* 2 dt*, je* 3 dt*] T (4.ie) 

and the weighting matrices \Q\ and [R] contain associated values presented in Sec. 4.1. The matrices 
[A], [Q], and [B] are dimensioned 9x9, 9x9, and 9x6, respectively, in accordance with the number 
of regulated state variables for this control scheme. Figure 4.19 shows the eigenvalues of the matrix 
[A] — [B] [ K ], or the closed-loop poles of the system. The diagonal lines mark the points at which 
a damping ratio of 1 /\/2 is obtained. The transpose of the 6x9 gain matrix [ K ] associated with 
these poles is given in Eq. (4.17). 
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Attitude-Hold Closed-Loop Poles 



Figure 4.19: Attitude- Hold Closed-Loop Poles 
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(4.17) 


The simulation results presented in Figs. 4.20-4.24 are obtained with this control law, together 
with the initial values of the state variables employed in Sec. 4.5. The commanded attitude for 
this simulation is = (9g = ^3 = 0; hence, a small- angle reorientation is necessary at first to 
align the spacecraft with the commanded attitude. Following this, the control law counteracts the 
disturbance torques to maintain that attitude until actuator saturation. 

A time history of the orientation of B in L, described with a body-three, 2-3-1 sequence, is 
shown in the upper plot of Fig. 4.20; the solid curve is used for 6 \ (pitch), the dashed curve for 
O 2 (yaw), and the dash-dot curve for 83 (roll). The lower plot of Fig. 4.20 illustrates the behavior 
of u\, v- 2 , and 113 , shown with solid, dotted, and dashed curves respectively. These results show 
the response due to the difference between the initial attitude and the commanded attitude. The 
attitude response reaches steady-state after roughly 1.5 orbits. 

Figure 4.21 displays the magnitude of momentum of the CMGs and flywheels in the upper and 
lower plots, respectively. Once again, the CMGs and flywheels are weighted evenly in the control 
law design, and as a result the momentum time histories are the same. The secular growth in the 
CMG and flywheel momentum is a consequence of the commanded attitude not being equal to the 
torque equilibrium attitude. 
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Figure 4.20: Attitude, Angular Velocity 
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Figure 4.21: Resultant Momentum Magnitude 
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Figure 4.24: Power, Power Error, Kinetic Energy Error 


Figure 4.22 shows the time histories of the CMG and flywheel control torques required to hold 
the attitude specified for this simulation. In each case, the torque component in the direction of 
bi is depicted with a solid curve, and those for b »2 and b 3 are indicated with dashed and dash-dot 
curves, respectively. It is apparent in this plot that the flywheel and CMG control torques are the 
same, which is to be expected because they are weighted evenly in the control law. In the steady- 
state, the control torques obtained with the TEA-seeking controller (see Fig. 4.5) are similar in 
magnitude to those observed here; however, a key difference is the presence of offsets or biases 
in the control torques for both types of actuator, caused by constant terms in the aerodynamic 
and gravity gradient disturbance torques. This offset exists for all three components of the control 
torque, leading to the secular momentum build up in the actuators; however it is most obvious in 
the b 2 direction. 

Figure 4.23 displays the angular speeds of the 6 flywheel rotors relative to B. Speeds of the 
pair whose spin axes are parallel to bi are contained in the upper plot, with 114 represented by 
the solid curve, and U 5 indicated by the dashed. Speeds of the pairs whose spin axes are parallel 
to b 2 and t >3 are shown in the middle and lower plots, respectively. The behavior illustrated in 
this plot is clearly not acceptable; the maximum speed is exceeded quickly as the flywheels store 
momentum. The rotor speeds in each pair gradually approach each other in order to produce the 
maximum amount of torque in one direction, which is needed to counter the secular external torque 
acting on the spacecraft. Just before the end of the third orbit, all flywheel pairs reach the same 
angular speeds, causing a singularity in the steering law that prevents the required power from 
being supplied, and thus the simulation is stopped at this point. Prior to this event, the power 
requirements and attitude control requirements are met simultaneously by the flywheel rotors, as 
is evident from Fig. 4.24. 
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Attitude-Hold with Flywheel Momentum Management Closed-Loop Poles 



Figure 4.25: Attitude-Hold Closed-Loop Poles, with Flywheel Momentum Management 


4.6.2 Attitude Hold with Flywheel Momentum Management 


In order to prevent the angular speeds of the flywheel pairs from becoming equal and causing a 
singularity in the steering law that renders the system incapable of providing the required power 
to the spacecraft, the control law may be reformulated by adding flywheel angular momenta and 
their integrals to the set of regulated state variables. Hence, 

{®} = [ 0 *, 02, 03, f Hf dt* , ■ ■ ■ , f 6\dt* , • • • ] T (4.18) 


In a sense, this configuration is a hybrid of the TEA-seeking and attitude-hold control laws, in 
which the CMGs are used to store momentum, and flywheels and CMGs both provide an attitude 
control torque. At the same time, the flywheel momentum feedback ensures that the angular speeds 
remain in the desired range of operation. The poles of the closed-loop system are shown in Fig. 4.25. 
The transpose of the 6 x 15 gain matrix [. K ] associated with these poles is given in Eq. (4.19). 
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-0.5997 

0.0000 

0.0000 

0.5007 

- 0.0000 

0.5133 

0.0000 

-0.2593 

-0.4069 

0.0000 

0.4765 

0.0000 

-3.4646 

- 0.0000 

0.0000 

-4.1502 

- 0.0000 

0.6967 

- 0.0000 

-0.4564 

0.5354 

- 0.0000 

-0.4644 

-0.3477 

- 0.0000 

-0.5374 

-0.7276 

- 0.0000 

-0.4741 


(4.19) 


55 


The attitude response produced with this control law is shown in Fig. 4.26, and is quite similar 
to the behavior obtained without flywheel momentum feedback, presented in Fig. 4.20. 

In contrast to the unacceptable motion exhibited in Fig. 4.23, the angular speeds of the flywheel 
pairs shown in Fig. 4.27 remain well separated and within the stated limits as a consequence of 
flywheel momentum feedback, thereby avoiding singularities while providing attitude control torque 
and power management simultaneously. 

Figure 4.28 shows the magnitude of CMG momentum in the upper plot, and the magnitude 
of flywheel momentum in the lower plot. The flywheel momentum remains small, apart from 
initial transients, due to the momentum feedback terms in the control law. The CMG momentum, 
however, grows without bound because the commanded attitude differs from the torque equilibrium 
attitude. The CMG momentum magnitude slope resulting from this control law configuration is 
essentially double that of Fig. 4.21 because the flywheels no longer aid the CMGs in storing the 
secular momentum. 

Figure 4.29 displays the control torque required to maintain the commanded attitude while 
simultaneously managing the flywheel momentum. Not surprisingly, the CMG control torque is 
larger than the flywheel control torque, and it is apparent upon comparing the plots that that 
CMGs are providing the torque required to null the effect of the constant terms in the disturbance 
torque, judging by the non-zero offsets. This offset results in the secular momentum buildup in 
the CMGs. The flywheels do not assist in cancelling the constant terms, but do aid the CMGs in 
overcoming the periodic terms in the disturbance resulting from aerodynamic and gravity torque. 
As in the previous example without flywheel momentum feedback, the offset is most apparent in 
the \)2 direction. 


Attitude of B in L 



Figure 4.26: Attitude, Angular Velocity 
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CMG Control Torque 



Figure 4.29: Control Torques 

4.6.3 Attitude Hold with Thrusters 

As shown in Secs. 4.6.1 and 4.6.2, holding a specified attitude other than the torque equilibrium 
attitude results in saturation of the momentum exchange devices. At that point some other means 
of controlling the attitude of the spacecraft is required; reaction control system (RCS) thrusters 
are typically used to provide the necessary attitude control torques. 

Attitude control thrusters are naturally discrete actuators, typically with a constant thrust 
magnitude and direction, and a variable thrust duration. The control laws developed for the 
CMGs and flywheels are continuous, therefore certain modifications must be made if thrusters are 
to be used. This implementation may be accomplished by specifying a thruster cycle time and 
holding the attitude control torque constant over this interval. The control system then determines 
the thruster on-time required to impart the same change in angular momentum to the spacecraft 
as the continuous control law demands. This objective may be accomplished optimally by choosing 
the thrusters in order to minimize the expenditure of propellant. As the thrusters operate with a 
certain minimum duration, and a maximum duration given by the cycle time, the discrete thruster 
torque may not always match the required continuous torque, typically for extremely small or 
extremely large torque commands. 

A minimum-fuel jet selection algorithm of this nature has been implemented using the thrust 
levels and directions, and jet locations given in Ref. [27], with a cycle time of 1 s. Since the 
flywheel system is only providing the power requirements, the flywheel angular speeds will remain 
in counter-rotation, as shown in Chap. 3, and an angular momentum feedback term is not required 
in this case. Therefore, the regulated state variables for this flight mode are given by Eq. (4.16), 
which results in the same state feedback gain matrix obtained in Sec. 4.6.1. 

Figure 4.30 shows the attitude response of this system, using the same initial conditions as the 
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previous examples, where the attitude control torques are provided by the thrusters rather than 
the CMGs or flywheels. The solid, dashed, and dash-dot lines represent the pitch, yaw and roll 
angles in the upper plot, and the generalized speeds u\, 112 , and 113 in the lower plot, respectively. 

The applied attitude control torque, tt, is shown in the upper plot of Fig. 4.31. Propellant is 
consumed at a rate of approximately 3.8 kg/orbit, as seen in the lower plot. 

4.6.4 Fuel Savings for Attitude-Hold Command 

Several control laws have been developed for using CMGs, flywheels, and thrusters to control the 
attitude of a spacecraft while simultaneously providing the required power. Holding a specified 
attitude other than the torque equilibrium attitude leads to saturation of the momentum exchange 
devices that provide the necessary attitude control torque without the benefit of momentum feed- 
back. RCS thrusters must then maintain the attitude until commanded to do otherwise (upon 
completion of a resupply docking, for example). Using the flywheels to aid the CMGs in controlling 
attitude may relieve some of the momentum burden which in turn will delay the saturation of these 
devices, thus reducing the propellant costs for holding a specified attitude for a certain period of 
time. Here we explore and quantify the potential fuel savings to be had in this manner. 

To quantify the propellant saving, the attitude of the spacecraft is held at 6 \ = 0* 2 = ^3 = 0 
for a period of one orbit, which is a typical scenario for rendezvous and docking procedures. In 
order to compare directly these results with material published in Ref. [28] regarding fuel costs for 
docking procedures, the inertia matrix of that report is used here. Specifically, 




' 94340560.0 

2316853.0 

5685949.0 ' 

jS/S*' 

= 1.36 x 

2316853.0 

79186536.0 

992227.0 



5685949.0 

992227.0 

147831360.0 


kg-m 2 


(4.20) 


The CMG saturation level is defined in Refs. [27] and [28] as 2,720 N-rn-s (2,000 ft-lbf-s). In this 
simulation the attitude angles and angular velocity of B relative to L all have initial values of zero, 
whereas those of the flywheel angular speeds and CMG momentum are the same as in previous 
examples. 

Figure 4.32 shows the fuel consumption resulting from application of the various control laws. 
The solid curve indicates the propellant expended when using the attitude-hold control law for 
which RCS thrusters are the sole actuators. This control method requires a propellant mass of 
4.21 kg to maintain the specified attitude for one orbit, which is comparable to the value of 4.7 kg 
obtained with a simulation of higher fidelity as reported on p. 17 of Ref. [28] for Stage 16A. 

The dashed line indicates the fuel consumption when the CMGs are used initially to provide 
the attitude control torques, until they saturate at 557 s, after which the RCS thrusters are used to 
hold the specified attitude for the remainder of the orbit. This requires 3.97 kg of propellant. Note 
that in this control configuration, the flywheels are only used to provide the required power and not 
to provide any control torques. Since the flywheels never saturate, they may be able to return the 
spacecraft to a TEA at the end of 1 orbit, and the CMGs can then desaturate themselves without 
the use of additional propellant. 

The dash-dot line displays the results of using the CMGs and flywheels to provide attitude 
control torques, with the benefit of flywheel momentum feedback. In this configuration, the CMG 
saturation is delayed slightly and occurs at 636 s. At that point the RCS thrusters are required 
to expend propellant in order to counteract the constant disturbance torques. The flywheels, 
however, are prevented from saturating by the momentum feedback portion of the control law, 
and are therefore available to aid the RCS thrusters in counteracting the periodic portion of the 
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Fuel Consumption During LVLH Hold 



Figure 4.32: Fuel Consumption 


disturbance torques, thus decreasing the propellant consumption even further. This method of 
maintaining the orientation of the spacecraft for one orbit results in an expenditure of 3.68 kg of 
fuel. The flywheels never saturate and therefore may be able to seek TEA at the end of 1 orbit, 
allowing the CMGs to desaturate themselves without the intervention of the RCS. 

The dotted line displays the result of removing the flywheel momentum feedback, so that the 
flywheels are used to store momentum along with the CMGs. This control mode results in a 
significant delay of CMG saturation until 813 s. Subsequently, the flywheels continue to store 
momentum until the maximum angular speeds of 60,000 rpm are reached after 1,997 s. The 
RCS thrusters are then used to provide the required attitude control torques. This configuration 
results in a significant amount of propellant savings, requiring only 2.06 kg of fuel; however, as the 
flywheels become saturated they can no longer store additional energy, and the power management 
task cannot be performed. This is clearly not an acceptable result unless the docking task can be 
completed prior to the flywheel saturation, in which case the maneuver can be completed without 
expending propellant. Because the CMG saturation level is set to 2,720 N-rn-s for this analysis, 
there is in fact a great deal of CMG momentum available at the end of 1 orbit; this may be enough 
to return the spacecraft to TEA and allow the CMGs to desaturate themselves without propellant 
expenditure. 

The results of the simulation performed to generate the dash-dot curve in Fig. 4.32 are shown in 
Figs. 4.33-4.37. Figure 4.33 illustrates the time histories of the attitude angles and inertial angular 
velocity of the spacecraft; disturbance torques lead the spacecraft on small excursions away from 
the initial, desired conditions. 

The flywheel angular speeds are seen in Fig. 4.34. It is clear that flywheel momentum feedback 
in the control law prevents rotor speed saturation and singularity in the steering law. 
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CMG Momentum Magnitude 



Figure 4.35: Resultant Momentum Magnitude 
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Figure 4.36: Control Torques 
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Figure 4.37: Power, Power Error, Kinetic Energy Error 


The magnitude of CMG and flywheel momentum is plotted in the upper and lower curves in 
Fig. 4.35, respectively. This plot displays the buildup of momentum and eventual saturation of the 
CMGs. Following saturation, the CMGs are no longer used to generate attitude control torques, and 
therefore the magnitude of the momentum remains constant. The CMG attitude control torques 
are then provided by the RCS thrusters, aided by those torques produced with the flywheel rotors. 

Figure 4.36 shows the time histories of the attitude control torques provided by the various 
actuators used in this simulation. The upper left plot contains the CMG control torque, which 
naturally goes to zero upon reaching the momentum saturation point. The RCS thrusters are then 
used to generate the required control torque, as shown in the upper right plot, which requires the 
consumption of propellant shown in the lower right plot. Flywheels are used to provide attitude 
control torque throughout the entire simulation; this torque is displayed in the lower left plot of 
Fig. 4.36. The flywheels are used also to generate the required power profile, which is tracked 
without error as shown in Fig. 4.37. 

The foregoing results demonstrate that in one orbit approximately 0.3 kg of propellant may be 
saved by using the flywheels to aid the CMGs in producing attitude control torques. Over 0.5 kg 
of propellant may be saved by using the combined CMG, flywheel, and RCS control law instead of 
relying solely upon the RCS thrusters. 

The simulation results shown here are for undamped flywheel rotors. Damping does inhibit 
the ability of the flywheels to provide the required power, as shown previously, but this effect is 
overcome by using the kinetic energy error feedback method. The damping does not, however, alter 
the fuel consumption results appreciably in the case we have studied. 
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Chapter 5 

Nonlinear Controllers 


The material of Chapter 4 deals with the control of attitude that is intended to be fixed, or nearly 
so, in a local-vertical-local-horizontal reference frame; it is therefore entirely reasonable to design 
and employ control laws based upon equations of motion linearized about such a condition, and 
expect them to perform well as long as attitude excursions remain small. The linear equations 
obtained previously provide a poor description of motion accompanying large-angle reorientations 
or slew maneuvers, hence, a nonlinear control design is required for such maneuvers. The purpose 
of the current chapter is to present first the design of an appropriate law, and then some results of 
its numerical exercise. 

Development of a globally asymptotically stable algorithm for reorientation is sketched out in 
Sec. 5.1. The approach adopted permits the imposition of state and control constraints similar to 
those obtained with the LQR technique by appropriate choices of the weighting matrices \Q\ and 
[R]. Section 5.2 contains results of a simulation in which the law directs CMGs to control attitude, 
while the only task of the flywheels is to furnish the required power. Consequences of allowing the 
CMGs to be assisted by the flywheels without and with the benefit of flywheel momentum feedback 
are reported in Secs. 5.3 and 5.4, respectively. 


5.1 Global Asymptotically Stable Reorientation 


The development and exercise of a controller for performing large-angle reorientation maneuvers is 
described in Refs. [29] and [30]. The controller is based upon feedback linearization and Lyapunov 
stability theory, and satisfies specified limits on slew rate and control torque. 

One of the motions that can be controlled by the algorithm is a change in the orientation of 
a rigid body B from one attitude fixed in any reference frame L, to another, where L moves in a 
Newtonian reference frame N. The angular velocities of B and L in N are related by Eq. (2.91), 


N u> B = N u L + L u B 


(5.1) 


and the dynamical equations governing the rotational motion of B in N can be expressed as 


B 


j B/B* _ 

dt 


L U B 


=«-( 


N 


u L + L u B 


fN 


) x • ( 


N 


u L + L u B ) 


_ J B/B* I _N, .L L 


(jJ — 


u B x N u> L 


(5.2) 


where \ B / B denotes the central inertia dyadic of B. The vector u represents the moment about 
B*, the mass center of B, of forces applied to B with the intent of controlling its rotational motion; 
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u is referred to simply as the control torque. With some manipulation, this relationship can be 
deduced from Eqs. (2.28) after setting p = 0 and letting u play the part of M, or one may refer to 
Eqs. (44) in Ref. [30]. It is worth noting that, unlike Eqs. (2.112)-(2.114) used to design the linear 
control laws, this relationship does not account for the presence of gravitational moment. 

During a large-angle reorientation it becomes possible for the second angle of an Euler angle 
sequence to take on a value of ±vr/2, leading to singular relationships between the time derivatives 
of the angles and angular velocity measure numbers, as can be seen in Appendix II of Ref. [16]. 
No such singularities occur in kinematical differential equations for quaternion elements, hence it 
is preferable to use them as kinematic variables instead of Euler angles, even though the later serve 
adequately for the linear controllers discussed in Chap. 4. 

Euler’s theorem on rotation states that every change in relative orientation of two bodies or 
reference frames L and B can be produced by means of a simple rotation of B in L. Quaternion 
elements qi , q 2 , <73, and (74 are defined in terms of a unit vector A fixed in L and B throughout a 
simple rotation having an angular displacement of 9. Referring to Eqs. (1.3.1)-(1.3.3) of Ref. [16], 

A . . 0 A .. 6 A . . 0 a 9 

<7i = Ai sin - , q 2 = A2SU1-, <73 = A 3 sm — , <74 = cos - 

where Ai, A2, and A3 are the dot products 

Aj = A • ij = A • b, (* = 1,2,3) (5.4) 

of A with li, I2, I3 and bi, b2, b3, dextral sets of orthogonal unit vectors fixed in L and B 
respectively, with 1* = bj (i = 1,2,3) prior to the rotation. The quaternion elements are not 
independent of each other, for they must obey the relationship 

<7i 2 + <72 2 + <73 2 + <74 2 = 1 (5.5) 


We have occasion in what follows to make use of the Euler vector q defined as 


A . . 9 
q = A sm - 

H 2 


(5.6) 


Kinematical differential equations for the quaternion can be expressed as in Eqs. (1.13.2) and 
(1.13.3) of Ref. [16], 


B 


4q = l(<MA, B + qx A, fl ) 

• 1 L B 

<74 — — 2 u • *4 

The control torque u derived in Refs. [29] and [30] takes the form 

n=( N u: L + L u B )x 1 B / B * ■ ( N u; L + L u; B )+l B / B * • V - L u B x 


(5.7) 

(5.8) 


- l B/B * • [A’(t)q + c( i ^ B ,q,q 4 ,t) L u B ] 


(5.9) 


where the terms on the first line of the right hand member are required for feedback linearization, 
and the second line represents feedback of the quaternion and the angular velocity of B relative to 
L. Thus, the dynamical equations of motion are reduced to the form 


B 


77 L “ B = -k(t)q - c ( L u B , q, <74, t) L uj B 


dt 


(5.10) 
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Global asymptotic stability is guaranteed for the dynamical system described by Eqs. (5.2), 
(5.7), and (5.8), and subject to the control torque u given by Eq. (5.9), when the scalar gain 
functions k(t) and c ( l uj b , q, < 74 , t) satisfy a set of conditions set forth in Ref. [29]. By asymptotic 
stability, we mean the magnitudes of L u> B and q are always driven to zero, and the scalar 54 always 
approaches ± 1 . Furthermore, satisfaction of a slew rate limit Cl 

L u B ■ L u B < Cl 2 t> 0 (5.11) 

is guaranteed when k and c meet additional conditions identified in Ref. [29] . A control constraint 
u 


u • u < u 2 t > 0 (5.12) 

is also satisfied when explicit gain functions k and c are constructed according to the specifications 
given in Ref. [30]. 

As indicated in Ref. [30] , design of the control law requires the analyst to specify the parameters 
c, the minimum value of the gain c, and k, the maximum value of gain k. A heuristic method is 
proposed to override the slew rate limit Cl with a smaller value in order to ensure that the magnitude 
of u remains less than or equal to u. This is done in Ref. [30] with the aid of two additional 
parameters, P 2 and ps, where a value of 0.5 is recommended for each. A nonlinear controller is 
constructed with the numerical values of these parameters reported in Table 5.1. 


Table 5.1: Nonlinear Controller Parameters 


Parameter 

Value 

c 

10~ 8 s ” 1 

k 

1 S ' 2 

P2 

0.5 

P3 

0.5 


The remainder of this chapter contains simulation results obtained by exercising this control 
law, where the role of the generic L in the foregoing development is played by the local-vertical- 
local horizontal reference frame described in Sec. 2.4. In all cases the initial attitude of B in L 
is given by 0\(to) = 50° (pitch), ^ 2 (^ 0 ) = 50° (yaw), and 0s(to) = 50° (roll), and the attitude 
specified for the conclusion of the reorientation is $1 = #2 = $3 = 0°. The initial and final 
values of the angular velocity L u B are of course 0 . Rotor spin speeds at the beginning are 
u^to) = ue(to) = us(to) = —20,000 rpm, and us(to) = = 'Ug(fo) = 20,000 rpm. Initial 

values of CMG momentum measure numbers are /ii(to) = ^ 2 (^ 0 ) = ^ 3 (^ 0 ) = 0- As before, changes 
in the mass distribution of the spacecraft are ignored; central principal moments of inertia are 
obtained from the moments and products given in Eq. 4.3, and bi, b 2 , and b 3 are taken to be 
parallel to central principal axes of inertia of S. The history of aerodynamic torque prescribed 
by Eq. 4.4 is appropriate when the attitude of this ISS configuration remains in the vicinity of 
TEA, but not when slew maneuvers are performed; therefore, no aerodynamic torque is applied in 
these simulations. As is the case in Chap. 4, gravitational moment is applied to the spacecraft in 
the simulation; since it is not included in Eq. (5.2), it presents an unmodeled disturbance to the 
controller. 
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Attitude of B in L 



1.5 


2.5 




Figure 5.1: Attitude 


5.2 Large-Angle Reorientation with CMGs 

Exploration of the performance of the nonlinear controller begins with attitude control performed 
solely by CMGs, while the flywheels provide the required power profile. Three values of the control 
torque limit u are considered: 1000 N-m, 500 N-m, and 350 N-m. A slew rate limit 0, specified 
as 0.35 deg/s, is automatically reduced by the aforementioned heuristic method for each value of 
u, yielding 0.085 deg/s, 0.052 deg/s, and 0.039 deg/s, respectively. The upper, middle, and lower 
plots of Figs. 5. 1-5.5 contain results associated with the high, middle, and low values of u. 

A time history of the orientation of B in L appears in Fig. 5.1; the solid curve is used for 6 \ 
(pitch), the dashed curve for 62 (yaw), and the dash-dot curve for 63 (roll). Simulations performed 
with values of u below 300 N-m give poor results (not shown) , because the controller is not given the 
authority to perform the feedback linearization. Naturally, as the control torque limit is decreased, 
the time to perform the maneuver is increased. 

Figure 5.2 shows the behavior of inertial angular velocity measure numbers u\, U 2 , and U 3 , 
presented with solid, dashed, and dash-dot curves respectively. 

Time histories of the CMG control torque corresponding to the three values of u are recorded in 
Fig. 5.3, where t±, 72 , and T 3 [having the same meanings as in Eqs. (2.93)— (2.95)] are displayed with 
solid, dashed, and dash-dot curves respectively. The associated control gains k and c are shown as 
functions of time in Fig. 5.4. This plot shows that the gains increase with time, until the maximum 
value of k is reached, and then the gains are held constant. Note that a decrease in the torque limit 
postpones the time at which the gains reach a maximum. 

the time at which the maximum gains are reached increases as the torque limit is decreased. 

Figure 5.5 displays the magnitudes of CMG momentum for the three levels of the control torque 
constraint. In each case it can be seen that the momentum capacity of the four ISS CMGs (18,981 
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Figure 5.2: Angular Velocity 
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Figure 5.3: Control Torques 
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Figure 5.6: Flywheel Angular Speeds 


N-m-s) is exceeded in less than one tenth of an orbit; as is well known, the CMGs are not adequate 
for performing large-angle reorientations or controlling the ISS at attitudes that differ significantly 
from a TEA. 


5.3 Large-Angle Reorientation with CMGs and Flywheels 

It stands to reason that matters can be improved, but only slightly, by allowing the flywheels and 
CMGs to share equally the burden of controlling attitude during a large-angle reorientation. This 
division of labor is accomplished by assigning half the nonlinear control torque u to each type of 
actuator, 


T = r = u/2 (5.13) 

The flywheels, of course, perform their duty of providing the required power as long as possible. 

Figure 5.6 shows the flywheel angular speeds associated with this control law configuration and 
u = 500 N-m. After slightly more than one tenth of an orbit, the rotor speeds in each flywheel pair 
become identical, causing a singularity in the steering law and preventing the required power from 
being supplied to the spacecraft. This state of affairs results from a lack of flywheel momentum 
feedback, as discussed in Sec. 4.6.1. 

5.4 Large-Angle Reorientation with CMGs, Flywheels, and Fly- 
wheel Momentum Feedback 

A simple way of making up for the lacuna in feedback is to augment the control with the product of 
flywheel momentum and integrals of momentum, and their corresponding LQR gains, an approach 
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similar to that taken in Sec. 4.6.2. 


where {x} = [ Hf , H. 
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> - [ K \ M 


(5.14) 


dt* ] T , and [K] is the LQR state feed- 
back gain matrix associated with {x}. Application of this control law, with u = 500 N-nr, to the 
reorientation maneuver gives rise to the results shown in Figs. 5.7-5.12. 

Time histories of the attitude and angular velocity of B in L are contained in Fig. 5.7, and the 
response is similar to what appears in the middle plots of Figs. 5.1 and 5.2. 

After inspection of the flywheel angular speeds plotted in Fig. 5.8, it is apparent that the 
steering law singularity has been avoided. In the steady state, the rotor speeds are similar to those 
illustrated in Fig. 4.27. 

The magnitudes of CMG and flywheel angular momenta are shown in Fig. 5.9. It is clear 
that the CMGs store the bulk of the momentum required to change the spacecraft’s orientation; 
unfortunately, the ISS CMGs would saturate less than one tenth of an orbit after the maneuver 
begins. After an initial period of transient behavior, flywheel momentum is kept near zero. 

Figures 5.10-5.11 show the control torques and nonlinear feedback gains, respectively. The 
CMG and flywheel control torques begin with the same initial values; however, their time histories 
are significantly different because the flywheel control torques are used to manage the flywheel 
angular momentum as well as assist the CMGs in providing the attitude control torques for the 
large-angle reorientation. The flywheels provide the required power profile, as demonstrated in 
Fig. 5.12. 

Augmenting the nonlinear control law as described in Eq. (5.14), and using the LQR gains 
obtained in Sec. 4.6.2, constitutes an ad-hoc approach at best; nevertheless, good performance is 
exhibited when conducting large-scale reorientations in L. The large amount of angular momentum 
required for such maneuvers, together with a capacity of current ISS CMGs that is limited even 
when supplemented with a significant flywheel energy storage system, dictates the use of RCS 
thrusters rather than momentum exchange devices. 

Stability properties of the augmented controller, formed by adding momentum feedback to the 
feedback linearization, have not been studied to date; this is a subject for future research. 
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Figure 5.9: Resultant Momentum Magnitude 
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Figure 5.10: Control Torques 
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Figure 5.11: Control Gains 
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Figure 5.12: Power, Power Error, and Kinetic Energy Error 







Chapter 6 


Summary of Results 


The purpose of this research is to obtain algorithms for simultaneous control of the orientation 
of an Earth-pointing spacecraft, the energy stored by counter-rotating flywheels, and the angular 
momentum of the flywheels and control moment gyroscopes (CMGs) used together as an integrated 
set of actuators for attitude control. Three important classes of motion relative to a local-vertical- 
local-horizontal reference frame L are examined. The first involves maintaining a torque equilibrium 
attitude (TEA), in which angular momentum of the actuators remains small and cyclic. In the 
second, attitude is held fixed, resulting in secular growth in angular momentum and eventual 
saturation of the momentum exchange devices. The third consists of a large-angle reorientation 
from one attitude, fixed in L, to another. 

In order to construct control laws and evaluate their performance with numerical simulations, 
dynamical equations of motion for a multibody spacecraft containing flywheels and CMGs are 
derived using Kane’s method, and expressed in vector-dyadic form. Each device is modeled as a 
rigid axisymmetric rotor; the spin axis of a flywheel is regarded as fixed in the spacecraft, whereas 
the spin axis of a CMG is not. The equations are nonlinear, and completely general with regard 
to the number and orientation of flywheel rotors, and the number of CMG rotors and gimbals. 
A set of twelve scalar equations is obtained by applying the generic relationships to the special 
case of a complex gyrostat consisting of a base body and three pairs of flywheels mounted in 
orthogonal directions. Existing literature contains equations for describing motion of a spacecraft 
with CMGs; they are shown to follow from the generic ones under two reasonable assumptions, 
namely that girnbal speeds are much less than the rotor spin speed, and reorientation of CMG 
rotors (and gimbals) does not significantly redistribute system mass. The exact equations for the 
complex gyrostat, and the approximate relationships associated with CMGs are combined to form 
approximate equations for a spacecraft with flywheels and CMGs, and subsequently linearized and 
nondimensionalized in preparation for design of linear control laws. 

We develop two steering laws whose purpose is to determine the flywheel motor torques neces- 
sary to meet the attitude control and power management requirements simultaneously. The first 
of these involves formation of a pseudo-inverse to solve an underdetermined system of equations by 
minimizing the sum of the squares of flywheel motor torques, whereas the second divides the power 
requirements evenly among the three pairs of rotors, resulting in a uniquely determined solution. 

Several laws are designed with the Linear Quadratic Regulator technique to control the first two 
types of motion described heretofore, seeking and maintaining TEA, and holding an attitude fixed. 
Controller performance is illustrated through numerical simulations involving the International 
Space Station (ISS). Energy storage is shown to be affected adversely by damping of the flywheel 
rotors, which must be expected under realistic conditions, but the problem is eliminated by feedback 
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of the error in flywheel rotational kinetic energy. A law that holds attitude without managing the 
angular momentum of flywheels or CMGs soon leads to undesirable flywheel rotor speeds and 
associated singularities in the steering laws, a situation that is rectified by the inclusion of flywheel 
momentum management. In order to continue holding attitude after the CMGs become saturated, 
we formulate a control law that uses reaction control system (RCS) thrusters and minimum-fuel jet 
selection logic. Thrusters, CMGs, and flywheels are used in various combinations to hold attitude 
for one orbit, and a significant reduction in propellant expenditure is possible when flywheels aid 
the CMGs. 

The third type of motion is addressed by the construction of nonlinear control laws for perform- 
ing large-angle reorientation maneuvers with CMGs and flywheels. An ad-hoc scheme for managing 
the flywheel angular momentum is proposed, similar in effect to what is obtained with linear control 
laws. This method performs well; however, proof of stability is a topic for further investigation. 
These large-angle maneuvers require a significant amount of momentum storage and as a result, 
RCS thrusters are required to provide the bulk of the necessary control torque. 

One topic of future research that may be valuable is the study of controller performance in 
the face of imperfect knowledge of spacecraft inertia properties, or flywheel rotor pairs whose axial 
moments of inertia are not identical. Inclusion of a pre-filter would allow the control laws to 
estimate these parameters from attitude response, and compensate for any significant changes in 
the mass distribution of the spacecraft over time. Another topic with possible merit is exploration 
of the best way to apportion the control torque between the flywheels and CMGs; an adaptive 
feature for distributing the workload unequally could be preferable to the practice of sharing it 
equally. Future spacecraft are likely to rely solely on flywheels, rather than a mixture of flywheels 
and CMGs, and the control laws developed here are easily applied to this special case. 

The promising results demonstrated here for control of power, momentum, and attitude of 
Earth-pointing spacecraft suggest that it would be worthwhile to examine the very important class 
of inertially oriented spacecraft. 
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